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Abstract. The Hamiltonian H specifies the energy levels and time evolution of 
a quantum theory. A standard axiom of quantum mechanics requires that H be 
Hermitian because Hermiticity guarantees that the energy spectrum is real and that 
time evolution is unitary (probability-preserving) . This paper describes an alternative 
formulation of quantum mechanics in which the mathematical axiom of Hermiticity 
(transpose + complex conjugate) is replaced by the physically transparent condition 
of space-time reflection (VT) symmetry. If H has an unbroken VT symmetry, then 
the spectrum is real. Examples of T^T-symmetric non-Hermitian quantum-mechanical 
Hamiltonians are H = + ix^ and H ^ — x^. Amazingly, the energy levels of these 
Hamiltonians are all real and positive! 

Does a PT-symmetric Hamiltonian H specify a physical quantum theory in which 
the norms of states are positive and time evolution is unitary? The answer is that if H 
has an unbroken VT symmetry, then it has another symmetry represented by a linear 
operator C. In terms of C, one can construct a time- independent inner product with 
a positive-definite norm. Thus, T'T-symmctric Hamiltonians describe a new class of 
complex quantum theories having positive probabilities and unitary time evolution. 

The Lee Model provides an excellent example of a T'T-symmetric Hamiltonian. 
The renormalizcd Lee-model Hamiltonian has a negative-norm "ghost" state because 
renormalization causes the Hamiltonian to become non-Hermitian. For the past 50 
years there have been many attempts to find a physical interpretation for the ghost, 
but all such attempts failed. The correct interpretation of the ghost is simply that the 
non-Hermitian Lee Model Hamiltonian is T'T-symmetric. The C operator for the Lee 
Model is calculated exactly and in closed form and the ghost is shown to be a physical 
state having a positive norm. The ideas of VT symmetry are illustrated by using many 
quantum-mechanical and quantum-field-theoretic models. 
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1. Introduction — New Kinds of Quantum Theories 

The theory of quantum mechanics is nearly one hundred years old and because there have 
been so many experimental verifications of its theoretical predictions, it has become an 
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accepted component of modern science. In an introductory course on quantum physics, 
one learns the fundamental axioms that define and characterize the theory. All but one 
of these axioms are physical requirements. For example, the energy spectrum is required 
to be real because all measurements of the energy of a system yield real results. Another 
axiom requires that the energy spectrum be bounded below so that the system has a 
stable lowest-energy state. Yet another axiom requires that the time evolution of a 
quantum system be unitary (probability-conserving) because the expected result of a 
probability measurement of a state cannot grow or decay in time. A quantum theory 
of elementary particles must also satisfy the physical axioms of Lorentz covariance and 
causality. However, there is one axiom that stands out because it is mathematical rather 
than physical in character, and this is the requirement that the Hamiltonian which 
is the operator that expresses the dynamics of the quantum system, be Hermitian. 

The requirement that H be Hermitian dates back to the early days of quantum 
mechanics. The Hermiticity of H is expressed by the equation 

H = H\ (1) 

where the Dirac Hermitian conjugation symbol f represents the combined operations of 
matrix transposition and complex conjugation. The mathematical symmetry condition 
dl]) is physically obscure but very convenient because it implies that the eigenvalues of 
H are real and that the time-evolution operator e~*^* is unitary. 

Hamiltonians that are non-Hermitian have traditionally been used to describe 
dissipative processes, such as the phenomenon of radioactive decay. However, these 
non-Hermitian Hamiltonians are only approximate, phenomenological descriptions of 
physical processes. They cannot be regarded as fundamental because they violate the 
requirement of unitarity. A non-Hermitian Hamiltonian whose purpose is to describe 
a particle that undergoes radioactive decay predicts that the probability of finding the 
particle gradually decreases in time. Of course, a particle cannot just disappear because 
this would violate the conservation of probability; rather, the particle transforms into 
other particles. Thus, a non-Hermitian Hamiltonian that describes radioactive decay 
can at best be a simplified, phenomenological, and nonfundamental description of the 
decay process because it ignores the precise nature of the decay products. In his book 
on quantum field theory Barton gives the standard reasons for why a non-Hermitian 
Hamiltonian cannot provide a fundamental description of nature [1] : "A non-Hermitian 
Hamiltonian is unacceptable partly because it may lead to complex energy eigenvalues, 
but chiefly because it implies a non-unitary S matrix, which fails to conserve probability 
and makes a hash of the physical interpretation." 

The purpose of this paper is to describe at an elementary level the breakthroughs 
that have been made in the past decade which show that while the symmetry condition 
([1]) is sufficient to guarantee that the energy spectrum is real and that time evolution 
is unitary, the condition of Dirac Hermiticity is not necessary. It is possible to describe 
natural processes by means of non-Hermitian Hamiltonians. We will show that the 
Hermiticity requirement ([1]) may be replaced by the analogous but physically transparent 
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condition of space-time reflection symmetry [VT symmetry) 

H = H'P'^ (2) 

witfiout violating any of the physical axioms of quantum mechanics. If H satisfies (^^, 
it is said to be VT symmetric. 

The notation used in this paper is as follows: The space-reflection operator, or parity 
operator, is represented by the symbol V. The effect of V on the quantum-mechanical 
coordinate operator x and the momentum operator p is to change their signs: 

VxV = —X and VpV = —p. (3) 

Note that "P is a linear operator and that it leaves invariant the fundamental 
commutation relation (the Heisenberg algebra) of quantum mechanics, 

xp — px = ih 1, (4) 

where 1 is the identity matrix. The time-reversal operator is represented by the symbol 
T. This operator leaves x invariant but changes the sign of p: 

TxT = X and TpT = —p. (5) 

Like the parity operator V, the time-reversal operator T leaves the commutation relation 
dl]) invariant, but this requires that T reverse the sign of the complex number i: 

TiT = -i. (6) 

Equation demonstrates that T is not a linear operator; T is said to be antilinear. 
Also, since V and T are reflection operators, their squares are the unit operator: 

V^ = r^ = 1. (7) 

Finally, the V and T operators commute: 

VT -TV = 0. (8) 

In terms of the V and T operators, we define the "PT-reflected Hamiltonian H^^ in ([2]) 
as H^^ = (VT)H{VT). Thus, if a Hamiltonian is VT symmetric [that is, if it satisfies 
([2])], then the VT operator commutes with H: 

H{VT) - {VT)H = 0. (9) 

A PT-symmetric Hamiltonian need not be Hermitian; that is, it need not satisfy 
the Hermiticity symmetry condition ([1]). Thus, it is possible to have a fully consistent 
quantum theory whose dynamics is described by a non-Hermitian Hamiltonian. Some 
examples of such non-Hermitian PT-symmetric Hamiltonians are 

H = p'^ + ix^, (10) 

and 

H = p'^-x\ (11) 

It is amazing indeed that the eigenvalues of these strange-looking Hamiltonians are all 
real and positive and that these two Hamiltonians specify a unitary time evolution even 
though they are non-Hermitian. 
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The Hamiltonians in fllOl) and fill I) are special cases of the general parametric family 
of PT-symmetric Hamiltonians 

H = p^ + x^{ix)% (12) 

where the parameter e is real. These Hamiltonians are all VT symmetric because they 
satisfy the condition in ([2]). It was shown in 1998 that when e > all of the eigenvalues 
of these Hamiltonians are entirely real and positive, but when e < there are complex 
eigenvalues [2]. We say that e > is the parametric region of unbroken VT symmetry 
and that e < is the parametric region of broken PT symmetry. (See Fig. [TJ) 

One can think of the non-Hermitian Hamiltonians in f|T2l) as complex extensions of 
the harmonic oscillator Hamiltonian H = + x'^. Indeed, the quantum theories defined 
by H are complex extensions of the conventional quantum theory of the harmonic 
oscillator into the complex domain. The general constructive principle that we are 
using in f|T2|) is to start with a Hamiltonian that is both Hermitian and VT symmetric. 
One then introduces a real parameter e in such a way that as e increases from the 
Hamiltonian is no longer Hermitian but its VT symmetry is maintained. One need 
not start with the harmonic oscillator. One can, for example, begin with any of the 
Hermitian Hamiltonians H = + x^^, where = 1, 2, 3, . . ., and introduce the 




Figure 1. Energy levels of the Hamiltonian H = + x'^{ixY as a function of the 
real parameter e. There are three regions: When e > 0, the spectrum is real and 
positive and the energy levels rise with increasing e. The lower bound of this region, 
e = 0, corresponds to the harmonic oscillator, whose energy levels are i5„ — 2n + 1. 
When — 1 < e < 0, there are a finite number of real positive eigenvalues and an infinite 
number of complex conjugate pairs of eigenvalues. As e decreases from to —1, the 
number of real eigenvalues decreases; when e < —0.57793, the only real eigenvalue is 
the ground-state energy. As e approaches — 1+, the ground-state energy diverges. For 
e < —1 there are no real eigenvalues. When e > 0, the VT symmetry is unbroken, but 
when e < the VT symmetry is broken. 
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parameter e as follows: H = + x"^^ [ixy. [The Hamiltonian in (1121) is just the special 
case = 1.] The properties of these Hamiltonians are discussed in [3]. 

We emphasize that these new kinds of Hamiltonians define valid and consistent 
quantum theories in which the mathematical condition of Dirac Hermiticity H = 
has been replaced by the physical condition of VT symmetry, H = H^^ . The condition 
in ([2]) that the Hamiltonian is VT symmetric is a physical condition because V and T are 
elements of the homogeneous Lorentz group of spatial rotations and Lorentz boosts. The 
real Lorentz group consists of four disconnected parts |4J: (i) The first part, called the 
proper orthochronous Lorentz group, is a subgroup of the Lorentz group whose elements 
are continuously connected to the identity, (ii) The second part consists of all of the 
elements of the proper orthochronous Lorentz group multiplied by the parity operator 
V. (iii) The third part consists of all of the elements of the proper orthochronous Lorentz 
group multiplied by the time-reversal operator T. (iv) The fourth part consists of all of 
the elements of the proper orthochronous Lorentz group multiplied by VT. Note that 
parts (ii) - (iv) are not subgroups of the Lorentz group because they do not contain the 
identity element. These four parts of the Lorentz group are disconnected because there 
is no continuous path in group space from one part to another. 

When we say that Lorentz invariance is a physical requirement of a theory, what 
we really mean is that the theory must be invariant under Lorentz transformations 
belonging to the proper, orthochronous Lorentz group. We know that the physical 
world is not invariant under the full homogeneous Lorentz group because it has been 
demonstrated experimentally that there exist weak processes that do not respect parity 
symmetry and other weak processes that do not respect time-reversal symmetry. 

One can extend the real Lorentz group to the complex Lorentz group [4]. (To 
perform this extension it is necessary to make the crucial assumption that the eigenvalues 
of the Hamiltonian are real and bounded below.) The complex Lorentz group consists 
of two and not four disconnected parts. In the complex Lorentz group there exists a 
continuous path in group space from the elements of the real proper, orthochronous 
Lorentz group to the elements of part (iv) of the real Lorentz group. There also exists 
a continuous path in group space from the elements of part (ii) to the elements of 
part (iii) of the real Lorentz group. Thus, while we know that the world is not invariant 
under V reflection or under T reflection, we are proposing here to consider the possibility 
suggested by complex group analysis that a fundamental discrete symmetry of the world 
is VT symmetry, or space-time reflection symmetry. 

The most important consequence of the discovery that non-Hermitian VT- 
symmetric Hamiltonians can define acceptable theories of quantum mechanics is that 
we now can construct many new kinds of Hamiltonians that only a decade ago would 
have been rejected as being unphysical because they violate the axiom of Hermiticity. In 
this paper we will examine some of these new Hamiltonians and discuss the properties 
and possible physical implications of the theories defined by these Hamiltonians. So far, 
there have been no experiments that prove clearly and definitively that quantum systems 
defined by non-Hermitian PT-symmetric Hamiltonians do exist in nature. However, 



Making Sense of Non-Hermitian Hamiltonians 



6 



one should keep an open mind regarding the kinds of theories that one is wilhng to 
consider. Indeed, Gell-Mann's "totahtarian principle" states that among all possible 
physical theories "Everything which is not forbidden is compulsory." 

1.1. Presentation and Scope of this Paper 

Like many research areas in science, the study of non-Hermitian Hamiltonians having 
real spectra began in a haphazard and diffuse fashion. There are numerous early 
examples of isolated and disconnected discoveries of such non-Hermitian Hamiltonians. 
For example, in 1980 Caliceti et al, who were studying Borel summation of divergent 
perturbation series arising from classes of anharmonic oscillators, were astonished to find 
that the eigenvalues of an oscillator having an imaginary cubic selfinteraction term are 
real [5]. In the summer of 1993, when I was visiting GEN Saclay, I learned that Bessis 
and Zinn- Justin had noticed on the basis of numerical work that some of the eigenvalues 
of the cubic Hamiltonian in ( ITOl) seemed to be real, and they wondered if the spectrum 
might be entirely real [6]. [Their interest in the Hamiltonian (fTOl) was piqued by the 
Lee- Yang edge singularity.] In 1982 Andrianov, who was doing perturbative studies 
of —x^ potentials, found evidence that such theories might have real eigenvalues [7j. 
In 1992 Hollowood [8] and Scholtz et al. [9] discovered in their own areas of research 
surprising examples of non-Hermitian Hamiltonians having real spectra. 

The field of PT-symmetric quantum mechanics was established in 1998 with the 
discovery by Bender and Boettcher that the numerical conjectures of Bessis and Zinn- 
Justin were not only valid, but were just one instance of a huge class of non-Hermitian 
Schrodinger eigenvalue problems whose spectra are entirely real and positive [2] . Bender 
and Boettcher showed that the reality of the spectra was due to a symmetry principle, 
namely the condition of an unbroken space-time reflection symmetry, and they argued 
that this symmetry principle could replace the usual requirement of Dirac Hermiticity. 

This discovery by Bender and Boettcher relies on two essential mathematical 
ingredients. First, Bender and Boettcher used the techniques of analytic continuation 
of eigenvalue problems. These fundamental techniques were developed and used heavily 
in the early work of Bender and Wu on divergent perturbation series [TOl [TT] and were 
later used by Bender and Turbiner [12]. These techniques are crucial because they 
show how to analytically continue the boundary conditions of an eigenvalue problem 
as a function of a parameter in the Hamiltonian. Second, Bender and Boettcher used 
the delta-expansion techniques that had been discovered and developed by Bender et 
al. [I3] as a way to avoid divergent perturbation series. The delta expansion is a 
powerful perturbation-theory technique in which the small perturbation parameter is a 
measure of the nonlinearity of the problem. As a result, many crucial properties of the 
problem are exactly preserved as this parameter varies. In the case of PT-symmetric 
quantum mechanics, it is the reality of the eigenvalues that is exactly maintained as the 
perturbation parameter e in ( IT2l) is varied. 

Many researchers have contributed immensely to the development of PT-symmetric 
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quantum mechanics by discovering new examples and models, proving theorems, and 
performing numerical and asymptotic analysis. In the past few years a large and active 
research community has developed and there have been half a dozen international 
conferences on the subject of VT symmetry, pseudo-Hermiticity, and non-Hermitian 
Hamiltonians. The proceedings of these conferences provide a complete source of 
references [T^ [T5| [T6l ITTt IT8| [19] . By now, there have been so many contributions 
to the field that it is impossible to describe them all in this one paper. 

The purpose of this paper is to give an elementary introduction to this exciting and 
active field of research. In writing this paper, my hope is that the rate of new discoveries 
and the development of the field will continue at such a rapid pace that this review will 
soon become obsolete. 

1.2. Organization of this Paper 

In Sec. [2]we show that PT-symmetric Hamiltonians are complex extensions of Hermitian 
Hamiltonians and we discuss the key property of PT-symmetric Hamiltonians, namely, 
that their energy eigenvalues are real and bounded below. We discuss techniques for 
calculating eigenvalues. In Sec. [3] we discuss classical PT-symmetric Hamiltonians and 
show how to continue ordinary classical mechanics into the complex domain. The 
crucial theoretical questions of conservation of probability and unitary time evolution 
are addressed in Sec. HI In Sec. [5] we illustrate the theory of PT-symmetric quantum 
mechanics by using a simple 2x2 matrix Hamiltonian. Then in Sec. [6] we explain how to 
calculate the C operator, which is the central pillar of VT symmetry and which is needed 
to construct the Hilbert space for the theory. In the next two sections we discuss the 
many applications of PT-symmetric quantum theory. We discuss quantum-mechanical 
applications in Sec. [7] and quantum-field-theoretic applications in Sec. [HI Finally, in 
Sec. [9] we make some brief concluding remarks. 

2. Determining the Eigenvalues of a PT-Symmetric Hamiltonian 

In this section we show how to calculate the eigenvalues of a PT-symmetric Hamiltonian. 
We begin by pointing out that the Hamiltonian operator defines and determines the 
physical properties of a quantum theory in three important ways: 

(i) The Hamiltonian determines the energy levels of the quantum theory. To find these 
energy levels one must solve the time-independent Schrodinger eigenvalue problem 

= Eip. (13) 

This equation usually takes the form of a differential equation that must be solved 
subject to boundary conditions on the eigenfunction ip. In the case of a VT- 
symmetric Hamiltonian it is crucial that the boundary conditions be imposed 
properly. We emphasize that for a quantum theory to be physically acceptable 
the energy eigenvalues of the Hamiltonian must be real and bounded below. 
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(ii) The Hamiltonian specifies the time evolution of the states and operators of the 
quantum theory. To determine the time evolution of a state il){t) in the Schrodinger 
picture we must solve the time-dependent Schrodinger equation 

i^^m = nm- (14) 

The solution to this first-order differential equation is straightforward because H is 
assumed to be independent of time: 

m = e-^^V(O). (15) 

We call e"*''^* the time- evolution operator. In conventional quantum mechanics the 
time-evolution operator is unitary because the Hamiltonian H is Hermitian. As a 
result, the norm of the state ilj{t) remains constant in time. The constancy of the 
norm is an essential feature of a quantum system because the norm of a state is a 
probability, and this probability must remain constant in time. If this probability 
were to grow or decay in time, we would say that the theory violates unitarity. In 
PT-symmetric quantum mechanics H is not Dirac Hermitian, but the norms of 
states are still time independent. 

(iii) The Hamiltonian incorporates the symmetries of the theory. As an example, 
suppose that the Hamiltonian H commutes with the parity operator V. We then 
say that the Hamiltonian is parity invariant. Since "P is a linear operator, we 
know that the eigenstates of the Hamiltonian [the solutions to ( |T3i) ] will also be 
eigenstates of V. Thus, the eigenstates of H will have a definite parity; they will 
be either even or odd under space reflection. (Of course, a general state, which is 
a linear combination of the eigenstates of H, need not have definite parity.) 

2.1. Broken and Unbroken VT Symmetry 

For the case of a PT-symmetric Hamiltonian, the VT operator commutes with the 
Hamiltonian H [see (|9])]. However, VT symmetry is more subtle than parity symmetry 
because the VT operator is not linear. Because VT' is not linear, the eigenstates of H 
may or may not be eigenstates of VT' . 

Let us see what may go wrong if we assume that an eigenstate of the Hamiltonian 
H is also an eigenstate of the VT operator. Call the eigenvalue A and express the 
eigenvalue condition as 

VTi) = Xil). (16) 
Multiply ([I6D by VT on the left and use the property that {VTf = 1 [see (jTl) and (IHD]: 

^jJ = {VT)X{VTfilj. (17) 
Since T is antilinear [see ([6])], we get 

^ = =\\\^il,. (18) 

Thus, |Ap = 1 and the eigenvalue A of the VT operator is a pure phase: 

A = e^". (19) 
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Next, multiply the eigenvalue equation (IT^ by VT on the left and again use the 
property that {VTY = 1: 

{VT)H^ = {VT)E{Vrf^. (20) 

Using the eigenvalue equation (fT6l) and recalling that VT commutes with we get 

HXil) = {VT)E{Vr)Xij. (21) 
Finally, we again use the property that T is antilinear to obtain 

EXip = E*\^. (22) 

Since A is nonzero [see f|T9l) ]. we conclude that the eigenvalue E is real: E = E*. 

In general, this conclusion is false, as Fig. [1] clearly demonstrates. When e < 0, 
some of the eigenvalues have disappeared because they are complex. On the other 
hand, for the restricted region e > this conclusion is correct; all of the eigenvalues are 
indeed real. We are then led to make the following definition: If every eigenfunction 
of a PT-symmetric Hamiltonian is also an eigenfunction of the VT operator, we say 
that the VT symmetry of H is unbroken. Conversely, if some of the eigenfunctions of a 
PT-symmetric Hamiltonian are not simultaneously eigenfunctions of the VT operator, 
we say that the VT symmetry of H is broken. 

The correct way to interpret fl22l) is that if a Hamiltonian has an unbroken VT 
symmetry, then all of its eigenvalues are real. Thus, to establish that the eigenvalues of 
a particular PT-symmetric Hamiltonian are real, it is necessary to prove that the VT 
symmetry of H is unbroken. This is difficult to show, and it took several years after the 
discovery of the family of PT-symmetric Hamiltonians in (fT2l) before a complete and 
rigorous proof was finally constructed by Dorey et al. in 2001 ^20i,i21jl§| Many others have 
contributed to the rigorous mathematical development of the theory of VT symmetry. 
These include Shin [25], Pham [26], Delabaere [27], Trinh [28], Weigert [291 [31], and 
Scholtz and Geyer [32] . Mostafazadeh generalized VT symmetry to pseudo-Hermiticity 
(see Subsec. 14.51) . 

2.2. Boundary Conditions for the Schrddinger Eigenvalue Problem 

Our objective in this section is to show how to calculate the eigenvalues of the 
Schrodinger eigenvalue problem (fT3|) . The most direct approach is to write ( |T3|) as 
a differential equation in coordinate space. The principal conceptual difficulty that we 
face in solving this differential equation is in identifying and understanding the boundary 
conditions on the coordinate-space eigenfunctions. 

§ The proof by Dorey et al. draws from many areas of theoretical and mathematical physics and 
uses spectral determinants, the Bethe ansatz, the Baxter T-Q relation, the monodromy group, and 
an array of techniques used in conformal quantum field theory. This proof is too technical to be 
described in this paper, but it was a significant advance because it establishes a correspondence between 
ordinary differential equations and integrable models. This correspondence is known as the ODE-IM 
correspondence [HI [231 [M] • 
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To write (IT^ in coordinate space, we make the standard transcriptions 

X ^ X and p —> —i-r-, (23) 

ax 

except that we treat the variable x as complex. The Schrodinger eigenvalue problem 
f fT3|) then takes the form 

- i)"{x) + x^{ixyi){x) = Eipix). (24) 

Although we cannot solve this equation exactly for arbitrary e, we can easily find 
the possible asymptotic behaviors of its solutions by using the WKB approximation. In 
general, for any differential equation of the form —y"{x) + V{x)y{x) = 0, where V{x) 
is a function that grows as |a;| oo, we know that the exponential component of the 
asymptotic behavior of y{x) for large has the form 



y[x) ~ exp 



± J ds^/V{s) 



(25) 



To identify the appropriate boundary conditions to impose on '?/'(x), we consider first 
the Hermitian case e = (the harmonic oscillator). From fl25|) we can see immediately 
that the possible asymptotic behaviors of solutions are iIj{x) ~ exp (±ix^). The usual 
requirement that the eigenfunction be square-integrable implies that we must choose 
the negative sign in the exponent, and therefore the eigenf unctions are Gaussian-like 
for large This result extends into the complex- a; plane: If the eigenfunctions vanish 
exponentially on the real-x axis for large they must also vanish in two wedges of 
opening angle ^vr in the complex plane centered about the positive-real and negative-real 
axes. These wedges are called Stokes wedges [30] . 

What happens as e increases from 0? As soon as e > (e noninteger), a logarithmic 
branch point appears at the origin x = 0. Without loss of generality, we may choose 
the branch cut to run up the imaginary axis from x = to x = ioo. In this cut plane 
the solutions to the differential equation are single-valued. From the asymptotic 
behavior of ipi^) in (p5|) . we deduce that the Stokes wedges rotate downward into the 
complex-x plane and that the opening angles of the wedges decrease as e increases. 

There are many wedges in which tplx) —>■ as |x| —>■ oo. Thus, there are 
many eigenvalue problems associated with the differential equation ( !24l) . We choose 
to continue the eigenvalue differential equation (12^ smoothly away from the location 
of the harmonic oscillator wedges at e = 0. (A detailed description of how to extend 
eigenvalue equations into the complex plane may be found in Ref. [12] •) For e > the 
center lines of the left and right wedges lie at the angles 

evr e TT . . 

t^ieft = — vr H and o'neht = • (26) 

e + 4 2 ^ e + 4 2 ^ ' 

The opening angle of each of these wedges is A = 27r/(e + 4). The differential equation 

flM|) may be integrated along any path in the complex- a; plane as long as the ends of the 

path approach complex infinity inside the left wedge and the right wedge. Note that 

these wedges contain the real-x axis when — 1 < e < 2. However, as soon as e is larger 

than 2, the wedges rotate below the real-x axis. These wedges are shown in Fig. [2l 
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Figure 2. Stokes wedges in tlie complex- a; plane containing the contour on which 
the eigenvalue problem for the differential equation for e = 2.2 is posed. In these 
wedges 'ip{x) vanishes exponentially as |a;| — > oo. The eigenfunction il}{x) vanishes most 
rapidly at the centers of the wedges. 

Notice that the wedges in Fig. [2] are mirror images of one another if they are 
reflected through the imaginary- a; axis. This left-right symmetry is the coordinate- 
space realization of VT symmetry. If we choose any point x in the complex-x plane 
and perform a parity reflection, then x —x. Time reversal replaces i by —i as we 
saw in ([6]), and so T replace —x by its complex conjugate —x*. Thus, in the coordinate 
representation VT symmetry is left-right symmetry. 

2.3. The Flaw in Dyson's Argument 

The quantum theories that we are considering in this paper are obtained by extending 
real quantum mechanics into the complex domain, as explained in Subsec. 12. 2[ The 
notion of analytically continuing a Hamiltonian into the complex plane was first 
discussed in 1952 by Dyson, who argued heuristically that perturbation theory for 
quantum electrodynamics diverges [33]. Dyson's argument consists of rotating the 
electric charge e into the complex-e plane: e — ie. Applied to the standard quantum 
anharmonic-oscillator Hamiltonian H = + gx^, Dyson's argument goes as follows: 
Rotate the parameter g anticlockwise in the complex-^f plane from positive g to —g. Now, 
the potential term in the Hamiltonian is no longer bounded below, so the resulting theory 
has no ground state. Hence, the ground-state energy EQ^g) has an abrupt transition at 
(7 = 0, which implies that Eo{g) must have a singularity at g = 0. 

Following Dyson's reasoning, one would think that the spectrum of the Hamiltonian 
( ITTl) . which is obtained by setting e = 2 in (fT2l) . would not be bounded below, and one 
might conclude that this Hamiltonian is mathematically and physically unacceptable. 
However, this heuristic argument is flawed. While the ground-state energy of the 
quantum anharmonic oscillator does indeed have a singularity at g = 0, the spectrum 
of the Hamiltonian (fTTj) that is obtained by analytically continuing a parameter in 
the Hamiltonian remains ambiguous until the boundary conditions satisfied by the 
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eigenfunctions are specified. The term "bounded below" is inappropriate if it relies on 
an ordering relation applied to a complex potential because ordering relations cannot be 
used for complex numbers. The concern that the spectrum of H in ( ITTl) is not bounded 
below is unfounded because PT-symmetric boundary conditions on the Schrodinger 
equation prohibit the occurrence of negative eigenvalues. 

The eigenvalues of H in ffTTl) depend crucially on the history of how this negative- 
coupling-constant Hamiltonian constant is obtained. There are two different ways to 
obtain H in ( fTTj) : First, one can substitute g = \g\e^^ into H = + gx^ and rotate from 
9 = to 9 = TV. Under this rotation, the ground-state energy -Eo(fi') becomes complex. 
Clearly, Eo{g) is real and positive when g > and complex when (7 < oljj] Second, one 
can obtain ffTTj) as a limit of the Hamiltonian H in f|T2l) as e : 0-^2. The spectrum of 
this complex Hamiltonian is real, positive, and discrete, as is shown in Fig. [H 

How can the Hamiltonian ( ITTi) possess two such astonishingly different spectra? 
The answer lies in the boundary conditions satisfied by the eigenfunctions ip{x). In 
the first case, in which 9 = a.Tgg is rotated in the complex-^f plane from to vr, tplx) 
vanishes in the complex- a; plane as |x| —>■ 00 inside the wedges — 7r/3 < argx < and 
— 47r/3 < argx < — vr. These wedges are not VT -symmetric reflections of one another. 
In the second case, in which the exponent e ranges from to 2, iIj{x) vanishes in the 
complex-x plane as |x| — 00 inside the PT-symmetric pair of wedges — 7r/3 < argx < 
and — vr < arga; < — 27r/3. We emphasize that in this second case the boundary 
conditions hold in wedges that are symmetric with respect to the imaginary axis; these 
boundary conditions enforce the VT symmetry of H and are responsible for the reality 
of the energy spectrum. 

Illustrative example: The harmonic oscillator Hamiltonian 

H =p^ + up-x? (w > 0) (27) 

is an elementary model that illustrates the dependence of the eigenvalues on the boundary 
conditions imposed on the eigenfunctions. Let us assume that the eigenfunctions of H vanish 
on the real- a: axis as a; — > ±c». The eigenfunctions have the form of a Gaussian exp (— ^wx^) 
multiplied by a Hermite polynomial. Thus, the Stokes wedges are centered about the positive 
real-x axis and have angular opening ^tt. The eigenvalues i?„ are given exactly by the formula 

£;„ = (n+i)c^ (n==0, 1, 2, 3, ...). (28) 

Now suppose that the parameter w is rotated by 180° into the complex-w plane from the positive 
axis to the negative axis so that w is replaced by — w. This causes the Stokes wedges in the 
complex-x plane to rotate by 90° so that the eigenfunctions now vanish exponentially on the 
imaginary- a; axis rather than on the real- a; axis. Also, as a consequence of this rotation, the 
eigenvalues change sign: 

En = - (n + i) w (n = 0, 1, 2, 3, . . .). (29) 

Notice that under the rotation that replaces to by —uj the Hamiltonian remains invariant, and yet 
the signs of the eigenvalues are reversed! This shows that the eigenspectrum depends crucially on 
the boundary conditions that are imposed on the eigenfunctions. 

II Rotating from 6' = to 6* = — tt, we obtain the same Hamiltonian as in (fTTj) . but the spectrum is the 
complex conjugate of the spectrum obtained when we rotate from = to 9 — n. 
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Apart from the eigenvalues, there is yet another striking difference between the two 
theories corresponding to H in ffTTl) . The expectation value of the operator x in the 
ground-state eigenfunction ipQ^x) is given by 

(0|0) - j^dxi^lix) ' ^ ^ 

where C is a complex contour that lies in the asymptotic wedges described above. The 
value of (0|a:|0)/(0|0) for H in (ITT]) depends on the limiting process by which we obtain 
H . If we substitute g = g^e^^ into the Hamiltonian H = + gx^ and rotate g from 
^ = to 6' = TT, we find by an elementary symmetry argument that this expectation 
value vanishes for all g on the semicircle in the complex-;? plane. The expectation 
value vanishes because this rotation in the complex-^f plane preserves parity symmetry 
(x — > — x). However, if we define H in (ITT!) by using the Hamiltonian in (fT2l) and 
by allowing e to range from to 2, we find that this expectation value is nonzero. 
In fact, this expectation value is nonvanishing for all e > 0. On this alternate path 
VT symmetry (reflection about the imaginary axis, x —x*) is preserved, but parity 
symmetry is permanently broken. (We suggest in Subsec. 18.41 that, as a consequence of 
broken parity symmetry, one might be able to describe the dynamics of the Higgs sector 
by using a PT-symmetric —g(p^ quantum field theory.) 

2.4- Using WKB Phase- Integral Techniques to Calculate Eigenvalues 

Now that we have identified the boundary conditions to be imposed on the eigenfunctions 
of the PT-symmetric Hamiltonian ( fT2l) , we can use a variety of techniques to calculate 
the eigenvalues of this Hamiltonian. Not surprisingly, it is impossible to solve the 
differential-equation eigenvalue problem (|2^ analytically and in closed form except in 
two special cases, namely, for e = (the harmonic oscillator) and for e ^ oo (the 
PT-symmetric version of the square- well potential, whose solution is given in Ref. |34j). 
Thus, it is necessary to use approximate analytic or numerical methods. 

The simplest analytic approach uses WKB theory, which gives an excellent 
approximation to the eigenvalues when e > 0. The WKB calculation is interesting 
because it must be performed in the complex plane rather than on the real-x axis. The 
turning points x± are those roots of = x^{ixy that analytically continue off the real 
axis as e increases from 0. These turning points, 

x_=E^e'^^^-^), x+ = E^e-'^^^-^) , (31) 

lie in the lower-half (upper-half) x plane in Fig. [2] when e > (e < 0). 

The leading-order WKB phase-integral quantization condition is given by 

{n + l/2)7r = I ^ dx ^jE-x^iixy. (32) 

J X- 

When e > this path lies entirely in the lower-half x plane, and when e = (the case 
of the harmonic oscillator) the path lies on the real axis. However, when e < the path 
lies in the upper-half x plane and crosses the cut on the positive imaginary-x axis. In 
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Table 1. Comparison of the exact eigenvalues (obtained with Runge-Kutta) and the 



e n 


WKB result 

-^oxact 


in dSl). 

-EwKB e 
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-E'cxact 


-E'wKB 
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1.156 267072 


1.0943 2 





1.477 149 753 


1.3765 
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4.109 228 752 


4.0895 
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6.003 386 082 


5.9558 
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7.562 273 854 


7.5489 
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11.802 433 593 


11.7690 


3 


11.314421818 


11.3043 


3 


18.458 818 694 


18.4321 


4 


15.291553 748 


15.2832 


4 


25.791792 423 


25.7692 


5 


19.451529125 


19.4444 


5 


33.694 279 298 


33.6746 


6 


23.766 740 439 


23.7606 


6 


42.093 814 569 


42.0761 


7 


28.217524 934 


28.2120 


7 


50.937278 826 


50.9214 


8 


32.789 082 922 


32.7841 


8 


60.185 767651 


60.1696 


9 


37.469 824 697 


37.4653 


9 


69.795 703 031 


69.7884 



this case there is no continuous path joining the turning points. Hence, WKB fails when 
e < 0. 

When e > 0, we deform the phase-integral contour so that it follows the rays from 
a;_ to and from to 

•1 



(n+i) 

We then solve for En- 



2 sin 



TT 



e + 2 



E^' 



En-- 



r(| + ^)v^(^ + i) 



2<!+4 

e+4 



ds Vl - 



[n 



oo 



(33) 



(34) 



sin (^) r (1 + ^) 

This formula gives a very accurate approximation to the eigenvalues plotted in Fig. [1] 
and it shows, at least in the WKB approximation, that the energy eigenvalues of H in 
( |T2l) are real and positive (see Table [I]). We can, in addition, perform a higher-order 
WKB calculation by replacing the phase integral by a closed contour that encircles the 
path joining the turning points (see Refs. [3] and [30]). 

It is interesting that the spectrum of the real Ixl^"*"^ potential strongly resembles that 
of the x'^{ixy potential. The leading-order WKB quantization condition (accurate for 
e > —2) is like that in except that the factor of sin (^) is absent. However, 
as e oo, the spectrum of \x\'^~^'^ approaches that of the square- well potential 
[En = {n + l)^7r^/4], while the energies of the complex x'^{ixy potential diverge, as 
Fig. [1] indicates. The energies of the \x\^ potential are shown in Fig. [3l 



2.5. Numerical Calculation of Eigenvalues 

There are several highly accurate numerical techniques for computing the energy 
spectrum that is displayed in Fig. [H The simplest and most direct method is to integrate 
the Schrodinger differential equation ( !24l) using a Runga-Kutta approach. To do so, we 
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Figure 3. Energy levels of the Hamiltonian H = + \x\^ as a function of the real 
parameter P. This figure is similar to Fig. [U but the eigenvalues do not pinch off 
and go into the complex plane because the Hamiltonian is Hermitian. (The spectrum 
becomes dense at P = 0.) 



convert this complex differential equation to a system of coupled, real, second-order 
equations. The convergence is most rapid when we integrate along paths located at the 
centers of the Stokes wedges and follow these paths out to oo. We then patch the two 
solutions in each Stokes wedge together at the origin. This procedure, which is described 
in detail in Ref. [2], gives highly accurate numerical results. 

An alternative is to perform a variational calculation in which we determine the 
energy levels by finding the stationary points of the functional 



is a three-parameter class of PT-invariant trial wave functions [35]. The integration 

contour C used to define {H){a, b, c) must lie inside the wedges in the complex-x plane 
in which the wave function falls off exponentially at infinity (see Fig. [2]). Rather than 
having a local minimum, the functional has a saddle point in (a, b, c)-space. At this 
saddle point the numerical prediction for the ground-state energy is extremely accurate 
for a wide range of e. This method also determines approximate eigenf unctions and 
eigenvalues of the excited states of these non-Hermitian Hamiltonians. Handy used the 
numerical technique of solving the coupled moment problem [36]. This technique, which 
produces accurate results for the eigenvalues, is the exact quantum-mechanical analog 
of solving the Schwinger-Dyson equations in quantum field theory. 




(35) 



where 




(36) 
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2.6. The Remarkable Case of a VT -Symmetric —x^ Potential 

The PT-symmetric — Hamiltonian in ( fTTI) . which is obtained by setting e = 2 in ( fT2l) . 
is particularly interesting because it is possible to obtain the energy spectrum by using 
real analysis alone; that is, one can focus on real x only and avoid having to perform 
analysis in the complex-x plane. One way to proceed is to deform the integration 
contour shown in Fig. [2] to the upper edges of the wedges so that it lies entirely on 
the real-x axis. If this is done carefully, the exact eigenvalues for this potential can be 
obtained by solving the Schrodinger equation (p^ subject to the boundary conditions 
that the potential be reflectionless [37]. That is, an incoming incident wave from the 
left gives rise to an outgoing transmitted wave on the right, but no reflected wave on 
the left. (This observation may have consequences in cosmological models, as explained 
in Subsec. 18.71 ) 

Another way to proceed is to show that the eigenvalues of the non-Hermitian — 
Hamiltonian are identical with the eigenvalues of a conventional Hermitian Hamiltonian 
having a positive x^ potential. A number of authors have observed and discussed this 
equivalence [3 [381 EH]- Here, we use elementary differential-equation methods [30] to 
prove that the spectrum of the non-Hermitian PT-symmetric Hamiltonian 

H = ^P''-gx^ {g>Q) (37) 
is identical to the spectrum of the Hermitian Hamiltonian 

H = —p^ + 4gx^ -hx r^x {g>0). (38) 
2m V fn 

(We have included the dimensional constants m, g, and h because they help to 
elucidate the physical significance of the spectral equivalence of these two very different 
Hamiltonians.) To show that H in fl37p and H in fl38p are equivalent, we examine the 
corresponding Schrodinger eigenvalue equations 

tlj"{x) - gx^tpix) = Etpix) (39) 



2m 
and 




'^^ x + Agx^ 1 $(x) = E^{x). (40) 



m 



We begin by moving the complex integration contour for the Schrodinger equation 
(!39|) to the real axis. To do so, we parameterize the integration contour using 



X = -2iL^/l + iy/L, (41) 

where 

L = \[hy{mg)Y^' (42) 

and y is a real parameter that ranges from — oo to cxd. A graph of the contour in (HT!) 
is shown in Fig. m The transformed differential equation then reads 
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Figure 4. Stokes wedges in the lower-half complex- a; plane for the Schrodinger 
equation arising from the Hamiltonian H in ([571) . The eigenfunctions of H decay 
exponentially as |a;| oo inside these wedges. Also shown is the contour in (|^T|) . 



Next, we perform a Fourier transform defined by 

/oo 
dye-^y^/^f{y). (44) 
-oo 

By this definition the Fourier transforms of a derivative and a product are given by 

f'{y) -^ipRp)/^ and yf{y) -^ihf'ip). (45) 
Thus, the transformed version of (HHl) reads 



lQgL''h^"{p) + ( — ^ + 32gL^h ) 0'(p) 



- 16gL' (^1 - i£^'kp) = (46) 
We expand and simphfy the differential equation in fl46p and get 

'2mL 

+ (l^ " 5z " ^^^^') ^^^^ = ^'^^^^^ ^^^^ 

Next, we eliminate the one-derivative term in the differential equation (1471) to 
convert it to the form of a Schrodinger equation. To do so, we substitute 

0(p) = e^(f)<l>(p), (48) 



where 



F{p) = ^p~ TTT^r^—TT^P^- (49) 



h"^ 192gmL^h 
The resulting equation is 



- + I - ^] Hp) = eHp). (50) 



Last, we rescale fl50l) by substituting 



p = xL^y32mg (51) 
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and we obtain the Schrodinger differential equation fl4Up . The completes the proof 
and verifies that the eigenvalues of the two Hamiltonians fl37|) and fl38l) have identical 
eigenvalues. This demonstration of equivalence is exact; no approximations were made 
in this argument]^ 



2. 7. Parity Anomaly 

The proof in Subsec. 12.61 that the Hamiltonians in (1371) and (1381) are equivalent helps to 
clarify some of the physical content of the non-Hermitian "PT-symmetric Hamiltonian 
( 137|) . The interpretation of the result in ( HOl) is that the linear term in the potential of 
the equivalent quartic Hermitian Hamiltonian in (l38ll is a parity anomaly. In general, 
an anomaly is a purely quantum (nonclassical) effect that vanishes in the classical limit 
h ^ 0. There is no classical analog of an anomaly because Planck's constant h does not 
appear in classical mechanics. 

We refer to the linear term in ( l38il as a parity anomaly for the following reason: Even 
though the PT-symmetric Hamiltonian (!37|) is symmetric under the parity reflections 
defined in ([3]), if does not respect parity symmetry. The violation of parity symmetry 
is subtle because it is contained in the boundary conditions that the eigenfunctions 
of the associated Schrodinger equation must satisfy. Since these boundary conditions 
are given at |x| = oo, the violation of parity symmetry is not detectable in any finite 
domain in the complex-x plane. Classical motion is a local phenomenon. Therefore, a 
classical particle that is moving under the influence of this Hamiltonian (see Sec. [3]) will 
act as if it is subject to parity-symmetric forces; the classical particle cannot feel the 
influences of quantum boundary conditions imposed at |x| = oo. In contrast, a quantum 
wave function is inherently nonlocal because it must obey boundary conditions that are 
imposed at |x| = oo. Thus, only a quantum particle "knows" about the violation of 
parity symmetry. To establish the equivalence between the Hamiltonians in (!37j) and 
fl38|) it was necessary to perform a Fourier transform [see (j^6l) ]. This transformation 
maps the point at x = oo to the point at p = 0, and this explains the presence of the 
linear parity- violating term in the potential of H in (1381) . The violation of parity is 
now a visible local effect in the Hamiltonian H . However, this violation of parity is 
proportional to h and evaporates in the classical limit h ^ 

The Hamiltonian fl38|) is Hermitian in the usual Dirac sense and its energy spectrum 
is bounded below. This Hamiltonian is also PT-symmetric because at every stage in 
the sequence of differential-equation transformations in Subsec. 12. 6[ VT symmetry is 
preserved. However, the variable x that gives rise to the parity anomaly in (HOl) is not 
a coordinate variable. Its behavior is that of a momentum variable because x changes 
sign under time reversal. 

^ A simple and exact transformation like the one presented in this subsection for mapping a VT- 
symmetric non-Hermitian Hamiltonian to a Hermitian Hamiltonian has been found only for the isolated 
case e = 2 in (fT2|) . [A more complicated spectral equivalence exists for the special case e = 4 (see 

Ref. [201 ED).] 
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The violation of parity symmetry at the quantum level has important physical 
implications. It is the lack of parity symmetry that implies that the one-point Green's 
function in the corresponding quantum field theory does not vanish. A possible 
consequence of this is that the elusive Higgs particle, which is a fundamental ingredient 
in the standard model of particle physics, is a quantum anomaly. In the next subsection, 
we show that the parity anomaly has a major impact on the spectrum of bound states 
in a quantum theory. 

2. 8. Physical Consequence of the Parity Anomaly: Appearance of Bound States in a 
VT -Symmetric Quartic Potential 

A direct physical consequence of the parity anomaly is the appearance of bound states. 
To elucidate the connection between the parity anomaly and bound states, we generalize 
the Hamiltonian ( 1371) to include a harmonic (x^) term in the potential: 

H = + - gx\ (52) 

The same differential-equation analysis used in Subsec. 12.61 straightforwardly yields the 
equivalent Hermitian Hamiltonian [38l [39] 

H = ^-h\[^x + Ag(x'-p^ . (53) 
2m y m \ 8g / 

Note that for these more general Hamiltonians the linear anomaly term remains 

unchanged from that in fl38l) . 

It was shown in an earlier paper ^\ that the Hamiltonian (|52|) exhibits bound 

states. In particle physics a bound state is a state having a negative binding energy. 

Bound states in the context of a quantum mechanics are defined as follows: Let the 

energy levels of the Hamiltonian be En {n = 0, 1, 2, . . .). The renormalized mass is the 

mass gap; that is, M = Ei — Eq. Higher excitations must be measured relative to the 

vacuum energy: En — Eq {n = 2, 3, 4, . . .). We say that the nth higher excitation is a 

bound state if the binding energy 

Bn = En-Eo- uM (54) 

is negative. If Bn is positive, then we regard the state as unbound because this state 
can decay into n 1-particle states of mass M in the presence of an external field. 

In Ref. [4lJ it was shown numerically that for small positive values of g the first 
few states of H in f l52|) are bound. As g increases, the number of bound states decreases 
until, when g/^^ is larger than the critical value 0.0465, there are no bound states at 
allEI 

Because H in (l52l) has the same spectrum as the Hermitian Hamiltonian in (1531) . it 
is easy to explain the appearance of bound states and to show that the bound states are 

+ In Ref. [41j a heuristic argument was given to explain why there is such a critical value. This 
argument is heuristic because the non-Hermitian Hamiltonian is evaluated for x in the complex plane. 
As is explained in Subsec. 12. 3( when x is complex, one cannot use order relationships such as > or <, 
which only apply to real numbers. 
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Figure 5. Potential of the Hermitian Hamiltonian ([55|) plotted as a function of the 
real variable x for the case e = and g = 0.046. The energy levels are indicated by 
horizontal lines. Because e = 0, there is no anomaly and the double-well potential is 
symmetric. Therefore, the mass gap is very small and thus there are no bound states. 

a direct consequence of the linear anomaly term. To probe the influence of the anomaly, 
we generalize (l53ll by inserting a dimensionless parameter e that measures the strength 
of the anomaly term: 

H=\p''-t^gx + Ag(^x^-^^ , (55) 

where for simplicity we have set m = fi = h = 1. 

If we set e = 0, there is no anomaly term and the potential is a symmetric double 
well. The mass gap for a double well is exponentially small because it is a result of the 
tunneling between the wells and thus the renormalized mass M is very small. Therefore, 
Bn in ( |5^ is positive and there are no bound states. In Fig. [5] we display the double- well 
potential and the first several states of the system for the case g = 0.046 and e = 0. 
There is a very small splitting between the lowest two states. 

If e = 1, the double- well potential is asymmetric and the two lowest states are not 
approximately degenerate. Therefore, bound states can occur near the bottom of the 
potential well. Higher-energy states eventually become unbound because, as we know 
from WKB theory, in a quartic well the nth. energy level grows like n'^^^ for large n. 
As g becomes large, the number of bound states becomes smaller because the depth of 
the double well decreases. For large enough g there are no bound states. In Fig. [6] we 
display the potential for e = 1 and for g = 0.046; for this g there is one bound state. 

Another way to display the bound states is to plot the value of the binding energy 
Bn as a function of n. For example, in Fig. [7] we display the bound states for e = 1 and 
g = 0.008333. Note that for these values there are 23 bound states. Observe also that 
the binding energy Bn is a smooth function of n. 

It is noteworthy that the bound-state spectrum depends so sensitively on the 
strength of the anomaly term in the Hamiltonian fl55l) . If e is slightly less than 1, 
the first few states become unbound, as shown in Fig. [HI In this figure g = 0.008333 and 
e = 0.9. If e is slightly greater than 1, the binding energy Bn is not a smooth function of 
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n for small n. In Fig. [Hlwe plot Bn as a function of n for g = 0.008333 and e = 1.1. Note 
that for these values of the parameters there are 30 bound states. Figures U\ El and [9] 
are strikingly different, which demonstrates the extreme sensitivity of the bound-state 
spectrum to the anomaly term. 

3. PT-Symmetric Classical Mechanics — The Strange Dynamics of a 
Classical Particle Subject to Complex Forces 

In this section we describe the properties of the PT-symmetric classical-mechanical 
theory that underlies the quantum-mechanical theory described by the Hamiltonian 
( IT2l) . We describe the motion of a particle that feels complex forces and responds by 
moving about in the complex plane. Several papers have been published in this area 
[HI |42l HH HH HHl |46] and we summarize here some of the surprising discoveries. 

One objective here is to explain heuristically how an upside-down potential like 
that in (fTTj) can have positive-energy quantum-mechanical eigenstates. One might think 
(incorrectly!) that since a classical particle would slide down the positive real axis 
to infinity, the corresponding quantum states would be unstable and the spectrum of 
the quantum system would be unbounded below. In fact, when e > (the region of 
unbroken VT symmetry), all but a set of measure zero of the possible classical paths are 
confined and periodic, and thus the classical particle does not slide off to infinity. When 
e < (the region of broken VT symmetry), the classical trajectories do indeed run off 
to infinity, and we can begin to understand why the energy levels of the corresponding 
quantum system are complex. 

The equation of motion of a classical particle described by H in ( 1121) follows from 
Hamilton's equations: 

dx dH ^ dp dH ... 
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Figure 6. Asymmetric potential well plotted as a function of the real variable x for the 
Hermitian Hamiltonian ([55]) with e = 1 and g — 0.046. The energy levels are indicated 
by horizontal lines. There is one bound state. The occurrence of bound states is due 
to the anomaly. 
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Figure 7. Binding energies Bn = En — Eq — nM plotted as a function of n for 
g = 0.008333 and € — \. A negative value of i?„ indicates a bound state. Observe 
that there are 23 bound states for these parameter values. Note that Bn is a smooth 
function of n. 
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Combining these two equations gives 

'^ = 2^{2 + e){^xf^\ (57) 

which is the complex version of Newton's second law, F = ma. 
We can integrate (157|) to give 

= ±VE + {^xr+^ (58) 

where E is the energy of the classical particle (the time- independent value of H). We 
treat time t as a real variable that parameterizes the complex path x{t) of this particle. 
Equation ( l58l) is a complex generalization of the concept that the velocity is the time 
derivative of the position (t> = ^). Here, t is real, but v and x are complex. 

We now describe and classify the solutions to Eq. (l58l) . Because the corresponding 
quantum theory possesses VT invariance, we restrict our attention to real values of E. 
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Figure 8. Binding energies Bn plotted as a function of n for (7 = 0.008333 and e = 0.9. 
The first five states have now become unbound and Bn is not a smooth function of n 
for n < 6. The next 12 states are bound, and in this region Bn is a smooth function 
of n. Comparison of this figure with Fig. [7] shows that the bound-state spectrum is 
exquisitely sensitive to the strength of the linear anomaly term. 
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Figure 9. Binding energies Bn plotted as a function of n for g = 0.008333 and e = 1.1. 
Note that there are 30 bound states and that Bn is not a smooth function of n when 
n is smah. 



Given this restriction, we can rescale x and t by real numbers so that without loss of 
generality Eq. (ISSjl reduces to 



3.1. The Case e = 

The classical solutions to Eq. ( j59l) have elaborate topologies, so we begin by considering 
some special values of e. For the simplest case, e = 0, there are two turning points and 
these lie on the real axis at ±1. To solve Eq. fl59l) we must specify the initial condition 
x(0). An obvious choice for x(0) is a turning point. If the path begins at ±1, there 
is a unique trajectory in the complex-x plane that solves (!59l) . This trajectory lies on 
the real axis and oscillates between the turning points. This is the usual sinusoidal 
harmonic motion. 

Choosing the energy determines the locations of the turning points, and choosing 
the initial position of the particle determines the initial velocity (up to a plus or minus 
sign) as well. So if the path of the particle begins anywhere on the real axis between the 
turning points, the initial velocity is fixed up to a sign and the trajectory of the particle 
still oscillates between the turning points. 

In conventional classical mechanics the only possible initial positions for the particle 
are on the real-x axis between the turning points because the velocity is real; all other 
points on the real axis belong to the so-called classically forbidden region. However, 
because we are analytically continuing classical mechanics into the complex plane, we can 
choose any point x(0) in the complex plane as an initial position. For all complex initial 
positions outside of the conventional classically allowed region the classical trajectory 
is an ellipse whose foci are the turning points. The ellipses are nested because no 
trajectories may cross. These ellipses are shown in Fig. [TOl 

The exact harmonic-oscillator (e = 0) solution to (l59l) is 




(59) 



x{t) = cos[arccosa;(0) ± t], 



(60) 
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where the sign of t determines the direction (clockwise or anticlockwise) in which the 
particle traces out the ellipse. For any ellipse the period is 27r. The period is the same 
for all trajectories because we can join the square-root branch points by a single finite 
branch cut lying along the real axis from x = —1 tox = 1. The complex path integral 
that determines the period can then be shrunk (by Cauchy's theorem) to the usual real 
integral joining the turning points. 

Note that all of the elliptical orbits in Fig. [10] are symmetric with respect to parity 
V (reflections through the origin) and time reversal T (reflections about the real axis) as 
well as to VT (reflections about the imaginary axis). Furthermore, V and T individually 
preserve the directions in which the ellipses are traversed. 

3.2. The Case e = 1 

When e = 1, there are three turning points. These turning points solve the equation 
ix^ = 1. Two lie below the real axis and are symmetric with respect to the imaginary 
axis: 

x_ = e"^^^/^ and x+ = e'^^/^ (61) 

Under VT reflection x_ and a;+ are interchanged. The third turning point lies on the 
imaginary axis at a;o = ^. 

Like the case e = 0, the trajectory of a particle that begins at the turning point x_ 
follows a path in the complex-x plane to the turning point at x+. Then, the particle 
retraces its path back to the turning point at x_, and it continues to oscillate between 
these two turning points. This path is shown on Fig. [TTJ The period of this motion 
is 2v^37rr (I) /r (I). A particle beginning at the third turning point xq exhibits a 




energy is i? = 1. The trajectories are nested ellipses with foci located at the turning 
points at x — ±1. The real line segment (degenerate ellipse) connecting the turning 
points is the conventional real periodic classical solution to the harmonic oscillator. 
All paths are closed orbits having the same period 2tt. 
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Figure 11. Classical trajectories in the complex-x plane for a particle of energy E = \ 
described by the Hamiltonian H = + ix^ . An oscillatory trajectory connects the 
turning points x±. This trajectory is enclosed by a set of closed, nested paths that 
fill the finite complex- a; plane except for points on the imaginary axis at or above the 
turning point = i- Trajectories that originate on the imaginary axis above x = i 
either move off to ioo or else approach xq, stop, turn around, and then move up the 
imaginary axis to ioo. 



completely distinct motion: It travels up the imaginary axis and reaches ioo in a finite 
time i/vFr (I) /r (I). This motion is not periodic. 

Paths originating from all other points in the finite complex-x plane follow closed 
periodic orbits. No two orbits may intersect; rather they are all nested, like the ellipses 
for the case e = 0. All of these orbits encircle the turning points x± and, by virtue 
of Cauchy's theorem, have the same period 2\/^T (|) /F (|) as the oscillatory path 
connecting x±. Because these orbits must avoid crossing the trajectory that runs up the 
positive imaginary axis from the turning point at xq = i, they are pinched in the region 
just below Xq , 8jS shown on Fig. [TTl 

3.3. The Case e = 2 

When e = 2, there are four turning points, two below the real axis and symmetric with 
respect to the imaginary axis, Xi = e~^*'^/^ and X2 = e~*'^/^, and two more located 
above the real axis and symmetric with respect to the imaginary axis, xs = e™!'^ and 
X4 = e^*'^/^. These turning points are solutions to the equation —x^ = 1. Classical 
trajectories that oscillate between the pair xi and X2 and the pair X3 and X4 are shown 
on Fig.[T2l The period of these oscillations is 2v^27rF (|) /F (|). Trajectories that begin 
elsewhere in the complex- a; plane are also shown on Fig. [TOl By virtue of Cauchy's 
theorem all these nested nonintersecting trajectories have the same period. All classical 
motion is periodic except for the special trajectories that begin on the real axis. A 
particle that begins on the real-x axis runs off to ±cxd; its trajectory is nonperiodic. 



Making Sense of Non-Hermitian Hamiltonians 



26 



2 



0.5 - 



1.5 - 



1 








-1.5 



-0.5 



-1 - 




-2 



-1.5 



-1 



-0.5 







0.5 



1.5 



Figure 12. Classical trajectories in the complex-x plane for a particle described by 
the Hamiltonian H = — x'^ and having energy E = 1. There are two oscillatory 
trajectories connecting the pairs of turning points xi and X2 in the lower-half x-plane 
and ^3 and 0:4 in the upper-half x-plane. [A trajectory joining any other pair of turning 
points is forbidden because it would violate VT (left-right) symmetry] The oscillatory 
trajectories are surrounded by closed orbits of the same period. In contrast to these 
periodic orbits, there is a special class of trajectories having unbounded path length 
and running along the real-x axis. 



3.4- Broken and Unbroken Classical VT Symmetry 

We can now understand heuristically why the energies of the corresponding VT- 
symmetric quantum systems are reaL In each of Figs. [THl [HI and [12] we can see that all 
of the orbits are localized and periodic. We may regard the pictured classical motions 
as representing particles confined to and orbiting around complex atomsl We can then 
use Bohr-Sommerfeld quantization to determine the discrete energies of the system: 



where C represents the orbit of a classical particle in the complex-x plane. By Cauchy's 
theorem, any closed orbit leads to the same result for the energy En, and because of 
VT symmetry the integral above gives a real value for the energy. 

The key difference between classical paths for e > and for e < is that in the 
former case the paths (except for isolated examples) are closed orbits and in the latter 
case the paths are open orbits. In Fig. [13] we consider the case e = —0.2 and display 
two paths that begin on the negative imaginary axis. Because e is noninteger, there is 
a branch cut and the classical particle travels on a Riemann surface rather than on a 
single sheet of the complex plane. In this figure one path evolves forward in time and 
the other path evolves backward in time. Each path spirals outward and eventually 
moves off to infinity. Note that the pair of paths forms a PT-symmetric structure. We 
remark that the paths do not cross because they are on different sheets of the Riemann 
surface. The function (zx)""'^ requires a branch cut, and we take this branch cut to lie 
along the positive imaginary axis. The forward-evolving path leaves the principal sheet 




(62) 
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Figure 13. Classical trajectories in the complex-x plane for the Hamiltonian in (|2T|) 
with e = —0.2. These trajectories begin on the negative imaginary axis very close to 
the origin. One trajectory evolves forward in time and the other goes backward in 
time. The trajectories are open orbits and show the particle spiraling off to infinity. 
The trajectories begin on the principal sheet of the Riemann surface; as they cross the 
branch cut on the positive imaginary axis, they visit the higher and lower sheets of the 
surface. The trajectories do not cross because they lie on different Riemann sheets. 



(sheet 0) of the Riemann surface and crosses the branch cut in the positive sense and 
continues on sheet 1. The reverse path crosses the branch cut in the negative sense and 
continues on sheet —1. Figure [13] shows the projection of the classical orbit onto the 
principal sheet. 

Figure [I3] shows why the energies of the quantum system in the broken VT- 
symmetric region e < are not real. The answer is simply that the trajectories are 
not closed orbits; the trajectories are open orbits, and all classical particles drift off to 
X = oo. As explained after (1321) . if we attempt to quantize the system for the case 
e < using the Bohr- Sommerf eld integral in (l62l) . the integral does not exist because 
the integration contour is not closed. 



3.5. Noninteger Values of e 

As e increases from 0, the turning points at x = 1 (and at shown in Fig. [10] 

rotate downward and clockwise (anticlockwise) into the complex-x plane. These turning 
points are solutions to the equation 1 + (ix)^"'"'^ = 0. When e is noninteger, this equation 
has many solutions that all lie on the unit circle and have the form 

/ 4]V-e\ 

X = exp ill— {N integer). (63) 



4 + 2e ^ 

These turning points occur in "PT-symmetric pairs (pairs that are symmetric when 
reflected through the imaginary axis) corresponding to the values (A^ = — 1, A^ = 0), 
{N = -2, N = 1), (A^ = -3, A^ = 2), (A^ = -4, A^ = 3), and so on. We label 
these pairs by the integer K {K = 0, 1, 2, 3, ...)so that the i^th pair corresponds to 
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Figure 14. Locations of tlie turning points for e = There are 11 T^T-symmetric 
pairs of turning points, with each pair being mirror images under reflection through 
the imaginary-x axis on the principal sheet. All 22 turning points lie on the unit 
circle on a five-sheeted Riemann surface, where the sheets are joined by cuts on the 
positive-imaginary axis. 

[N = —K — 1,N = K). The pair of turning points on the real-x axis for e = deforms 
continuously into the ii' = pair of turning points when e > 0. When e is rational, there 
are a finite number of turning points in the complex-x Riemann surface. For example, 
when e = y, there are 5 sheets in the Riemann surface and 11 pairs of turning points. 
The K = {] pair of turning points are labeled = — 1 and = 0, the K = 1 pair are 
labeled = —2 and = 1, and so on. The last {K = 10) pair of turning points are 
labeled A^ = — 11 and A^ = 10. These turning points are shown in Fig. [HI 

As e increases from 0, the elliptical complex trajectories in Fig. [10] for the harmonic 
oscillator begin to distort but the trajectories remain closed and periodic except for 
special singular trajectories that run off to complex infinity, as we see in Fig. [Til These 
singular trajectories only occur when e is an integer. Nearly all of the orbits that one 
finds are VT symmetric (left-right symmetric), and until very recently it was thought 
that all closed periodic orbits are VT symmetric. This is, in fact, not so [15]. Closed non- 
PT-symmetric orbits exist, and these orbits are crucial for understanding the behavior 
of the periods of the complex orbits as e varies. 
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Figure 15. A classical trajectory in the complex- cc plane for the complex Hamiltonian 
H = p'^x^{ixY~'^ . This compHcated trajectory begins at a;(0) = —T.li and visits 11 
sheets of the Riemann surface. Its period is approximately T = 255.3. This figure 
displays the projection of the trajectory onto the principal sheet of the Riemann 
surface. This trajectory does not cross itself. 




Figure 16. An enlargement of the classical trajectory x{t) in Fig. [15] showing the 
detail near the origin in the complex- a; plane. We emphasize that this classical path 
never crosses itself; the apparent self-intersections are paths that lie on different sheets 
of the Riemann surface. 



In Fig. [15] we display a PT-symmetric orbit of immense topological intricacy that 
visits many sheets of the Riemann surface. This figure shows a classical trajectory 
corresponding to e = vr — 2. The trajectory starts at x(0) = —7.1i and visits 11 sheets 
of the Riemann surface. Its period is T = 255.3. The structure of this orbit near the 
origin is so complicated that we provide a magnified version in Fig. [T6l 

The period of any classical orbit depends on the specific pairs of turning points that 
are enclosed by the orbit and on the number of times that the orbit encircles each pair. 
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As shown in Ref. |13], any given orbit can be deformed to a simpler orbit of exactly 
the same period. This simpler orbit connects two turning points and oscillates between 
them rather than encircling them. For the elementary case of orbits that enclose only 
the K = Q pair of turning points, the formula for the period of the closed orbit is 

The derivation of flMl) goes as follows: The period Tq is given by a closed contour integral 
along the trajectory in the complex-x plane. This trajectory encloses the square-root 
branch cut that joins the K = Q pair of turning points. This contour can be deformed 
into a pair of rays that run from one turning point to the origin and then from the 
origin to the other turning point. The integral along each ray is easily evaluated as a 
beta function, which is then written in terms of gamma functions. 

When the classical orbit encloses more than just the K = Q pair of turning points, 
the formula for the period of the orbit becomes more complicated [l3j. In general, there 
are contributions to the period integral from many enclosed pairs of turning points. We 
label each such pair by the integer j. The formula for the period of the topological class 
of classical orbits whose central orbit terminates on the i^'th pair of turning points is 




In this formula the cosines originate from the angular positions of the turning points in 
( l63l) . The coefficients aj{K, e) are all nonnegative integers. The jth coefficient is nonzero 
only if the classical path encloses the jth pair of turning points. Each coefficient is an 
even integer except for the j = K coefficient, which is an odd integer. The coefficients 
ttj {K, e) satisfy 

oo 

J2(^j{K,e) = k, (66) 

j=0 

where k is the number of times that the central classical path crosses the imaginary 
axis. Equation ( |66l) truncates the summation in (l65l) to a finite number of terms. 

The period Tq in flMl) of orbits connecting the K = turning points is a smoothly 
decreasing function of e. However, for classical orbits connecting the i^th {K > 0) pair 
of turning points, the classical orbits exhibit fine structure that is exquisitely sensitive 
to the value of e. Small variations in e can cause huge changes in the topology and in 
the periods of the closed orbits. Depending on e, there are orbits having short periods 
as well as orbits having long and possibly arbitrarily long periods. 

3. 6. Classical Orbits Having Spontaneously Broken VT Symmetry 

There is a general pattern that holds for all K. For classical orbits that oscillate between 
the Kih pair of turning points, there are three regions of e. The domain of Region I is 
< e < the domain of Region II is < e < AK, and the domain of Region III is 



Making Sense of Non-Hermitian Hamiltonians 



31 




Figure 17. Period of a classical trajectory beginning at the iV = 1 turning point in 
the complex-x plane. The period is plotted as a function of e. The period decreases 
smoothly for < e < 1 (Region I). However, when 1 < e < 4 (Region II), the period 
becomes a rapidly varying and noisy function of e. For e > 4 (Region III) the period 
is once again a smoothly decaying function of e. Region II contains short subintervals 
where the period is a small and smoothly varying function of e. At the edges of these 
subintervals the period suddenly becomes extremely long. Detailed numerical analysis 
shows that the edges of the subintervals lie at special rational values of e. Some of 
these special rational values of e are indicated by vertical line segments that cross the 
horizontal axis. At these rational values the orbit does not reach the N = —2 turning 
point and the VT symmetry of the classical orbit is spontaneously broken. 

AK < e. In Regions I and III the period is a small and smoothly decreasing function 
of e. However, in Region II the period is a rapidly varying and noisy function of e. We 
illustrate this behavior for the case K = 1 and K = 2 in Figs. [T7| and [181 

The abrupt changes in the topology and the periods of the orbits for e in Region II 
are associated with the appearance of orbits having spontaneously broken VT symmetry. 
In Region II there are short patches where the period is relatively small and is a slowly 
varying function of e. These patches are bounded by special values of e for which the 
period of the orbit suddenly becomes extremely long. Numerical studies of the orbits 
connecting the i^th pair of turning points indicate that these special values of e are 
always rational |16]. Furthermore, at these special rational values of e, the closed orbits 
are no^ PT-symmetric (left-right symmetric). Such orbits exhibit spontaneously broken 
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Figure 18. Period of a classical trajectory joining (except when VT symmetry is 
broken) the K = 2 pair of turning points. The period is plotted as a function of e. As 
in the K = I case shown in Fig. 1171 there are three regions. When < e < i (Re gion 
I), the period is a smooth decreasing function of e; when | < e < 8 (Region II), the 
period is a rapidly varying and choppy function of e; when 8 < e (Region III), the 
period is again a smooth and decreasing function of e. 



VT symmetry. Some special values of e at which spontaneously broken PT-symmetric 
orbits occur are indicated in Figs. [ITl and [18] by short vertical lines below the horizontal 
axis. These special values of e have the form where p is a multiple of 4 and q is odd. 

A broken- "PT-symmetric orbit is a failed "PT-symmetric orbit. Figure [T9l displays 
a spontaneously-broken-PT-symmetric orbit for e = |- The orbit starts at the N = 2 
turning point, but it never reaches the PT-symmetric turning point = —3. Rather, 
the orbit terminates when it runs into and is reflected back from the complex conjugate 
= 4 turning point [see fl63|) ]. The period of the orbit is short (T = 4.63). While this 
orbit is not VT (left-right) symmetric, it does possess complex-conjugate (up-down) 
symmetry. In general, for a non-PT-symmetric orbit to exist, it must join or encircle a 
pair of complex-conjugate turning points. 

If we change e slightly, VT symmetry is restored and one can only find orbits that 
are VT symmetric. For example, if we take e = 0.805, we obtain the complicated orbit 
in Fig. [20l The period of this orbit is large (T = 173.36). 

Broken-PT-symmetric orbits need not be simple looking, like the orbit shown in 
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Figure 19. A horseshoe-shaped non-T'T-symnietric orbit. This orbit is not symmetric 
with respect to the imaginary axis, but it is symmetric with respect to the real axis. 
The orbit terminates at a complex-conjugate pair of turning points. For this orbit 




Fig. [ini Indeed, they can have an elaborate topology. As an example we plot in 
Fig. [21] the complicated orbit that arises when e = ^. This orbit is a failed K = 
PT-symmetric orbit that originates at the = —4 turning point, but never reaches 
the PT-symmetric = 3 turning point. Instead, it is reflected back by the complex- 
conjugate A^ = —14 turning point. 

This study of classical orbits provides a heuristic explanation of the quantum 
transition from a broken to an unbroken VT symmetry as e increases past 0. The 
quantum transition corresponds to a change from open to closed classical orbits. 
Furthermore, we can now see why the quantum theory in the unbroken region is 
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Figure 21. Non-T'T-symmetric orbit for e = ^. This topologically complicated orbit 
originates at the TV = —4 turning point but does not reach the T^T-symmetric = 3 
turning point. Instead, it is reflected back at the complex-conjugate TV = — 14 turning 
point. The period of this orbit is T 186.14. 



unitary. At the classical level, particles are bound in a complex atom and cannot 
escape to infinity; at the quantum level the probability is conserved and does not leak 
away as time evolves. Quantum mechanics is obtained by summing over all possible 
classical trajectories, and in the case of "PT-symmetric classical mechanics we have seen 
some bizarre classical trajectories. To understand how summing over such trajectories 
produces PT-symmetric quantum mechanics will require much more research. 

4. PT-Symmetric Quantum Mechanics 

Establishing that the eigenvalues of many PT-symmetric Hamiltonians are real and 
positive raises an obvious question: Does a non-Hermitian Hamiltonian such as H in 
( IT2|) define a physical theory of quantum mechanics or is the reality and positivity of 
the spectrum merely an intriguing mathematical curiosity exhibited by some special 
classes of complex eigenvalue problems? Recall that a physical quantum theory must 
(i) have an energy spectrum that is bounded below; (ii) possess a Hilbert space of state 
vectors that is endowed with an inner product having a positive norm; (iii) have unitary 
time evolution. The simplest condition on the Hamiltonian H that guarantees that the 
quantum theory satisfies these three requirements is that H be real and symmetric. 
However, this condition is overly restrictive. One can allow H to be complex as long 
as it is Dirac Hermitian: H'^ = H. In this section we explain why we can replace the 
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condition of Hermiticity by the condition that H have an unbroken VT symmetry and 
still satisfy the above requirements for a physical quantum theoryEl 

4.1. Recipe for a Quantum- Mechanical Theory Defined by a Hermitian Hamiltonian 

For purposes of comparison, we summarize in this subsection the standard textbook 
procedure that one follows in analyzing a theory defined by a conventional Hermitian 
quantum-mechanical Hamiltonian. In the next subsection we repeat these procedures 
for a non-Hermitian Hamiltonian. 

(a) Eigenf unctions and eigenvalues of H. Given the Hamiltonian H one can write 
down the time-independent Schrodinger equation associated with H and calculate 
the eigenvalues and eigenf unctions ipni^)- Usually, this calculation is difficult 
to perform analytically, so it must be done numerically. 

(b) Orthogonality of eigenf unctions. Because H is Hermitian, the eigenfunctions of H 
will be orthogonal with respect to the standard Hermitian inner product: 



(7/^,0) = J dx[^{x)]*(f){x) (67) 

Orthogonality means that the inner product of two eigenfunctions ipm{x) and ipni^) 
associated with different eigenvalues 7^ E^ vanishes: 

i^m, (Pn) = 0. (68) 

(We do not discuss here the technical problems associated with degenerate spectra.) 

(c) Orthonormality of eigenfunctions. Since the Hamiltonian is Hermitian, the norm 
of any vector is guaranteed to be positive. This means that we can normalize the 
eigenfunctions of H so that the norm of every eigenfunction is unity: 

(7/^„,V^„) = 1- (69) 

(d) Completeness of eigenfunctions. It is a deep theorem of the theory of linear 
operators on Hilbert spaces that the eigenfunctions of a Hermitian Hamiltonian 
are complete. This means that any (finite-norm) vector x the Hilbert space can 

be expressed as a linear combination of the eigenfunctions of H: 

00 

X = X] an^Jn- (70) 

n=0 

The formal statement of completeness in coordinate space is the reconstruction of 
the unit operator (the delta function) as a sum over the eigenfunctions: 



00 

n=0 



[iPn{x)ri^n{y)=5{x-y). (71) 



* All of the ■pT-symmetric Hamiltonians considered in this paper are symmetric under matrix 
transposition. This matrix symmetry condition is not necessary, but it has the simplifying advantage 
that the we do not need to have a biorthogonal set of basis states. We can consider T'T-symmetric 
Hamiltonians that are not symmetric under matrix transposition, but only at the cost of introducing a 
biorthogonal basis [311 HZ] ■ 
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(e) Reconstruction of the Hamiltonian H and the Green's function G, and calculation 
of the spectral Zeta function. The Hamiltonian matrix in coordinate space has the 
form 

oo 

YM^n{x)ri;^{y)E^ = H{x,y) (72) 

n.=0 

and the Green's function is given by 

oo ^ 

YM'n{x)rMy)— = G{x,y). (73) 

The Green's function is the matrix inverse of the Hamiltonian in the sense that 

dyH{x,y)G{y,z) = 5{x-z). (74) 



The formula for the Green's function in (1731) allows us to calculate the sum of the 
reciprocals of the energy eigenvalues. We simply set x = y in fl73l) . integrate with 
respect to x, and use the normalization condition in fl69|) to obtain the result that 



oo ^ 

dxG{x,x) = '^^^—. (75) 



n=0 -"'^ 



The summation on the right side of f[75|) is called the spectral zeta function. This 
sum is convergent if the energy levels En rise faster than linearly with n. Thus, the 
spectral zeta function for the harmonic oscillator is divergent, but it exists for the 
Ixl^"*"^ and x'^{ixy potentials if e > 0. 

(f) Time evolution and unitarity. For a Hermitian Hamiltonian the time-evolution 
operator e~*^* [see ( fT5l) ] is unitary, and it automatically preserves the inner product: 

{x{t). Xit)) = {xiO)e'"\ e-^*x(t)) = (x(0), x(0)) . (76) 

(g) Observables. An observable is represented by a linear Hermitian operator. The 
outcome of a measurement is one of the real eigenvalues of this operator. 

(h) Miscellany. One can study a number of additional topics, such as the classical 
and semiclassical limits of the quantum theory, probability density and currents, 
perturbative and nonperturbative calculations, and so on. We do not address these 
issues in depth in this paper. 



4-2. Recipe for VT- Symmetric Quantum Mechanics 

Let us follow the recipe outlined in Subsec. 14.11 for the case of a non- Hermitian VT- 
symmetric Hamiltonian having an unbroken VT symmetry. For definiteness, we will 
imagine that the non- Hermitian Hamiltonian has the form in ( fT2|) . The novelty here 
is that we do not know a priori the definition of the inner product, as we do in the 
case of ordinary Hermitian quantum mechanics. We will have to discover the correct 
inner product in the course of our analysis. The inner product is determined by 
the Hamiltonian itself, so PT-symmetric quantum mechanics is a kind of "bootstrap" 
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theory. The Hamiltonian operator chooses its own Hilbert space (and associated inner 
product) in which it prefers to hve! 

(a) Eigenfunctions and eigenvalues of H. In Sec. [2] we discussed various techniques for 
determining the coordinate-space eigenfunctions and eigenvalues of a non-Hermitian 
Hamiltonian. We assume here that we have found the eigenvalues En by using either 
analytical or numerical methods and that these eigenvalues are all real. [This 
is equivalent to assuming that the VT symmetry of H is unbroken; that is, all 
eigenfunctions ipn{x) of H are also eigenfunctions of VT.]. 

(b) Orthogonality of eigenfunctions. To test the orthogonality of the eigenfunctions, we 
must specify an inner product. (A pair of vectors can be orthogonal with respect 
to one inner product and not orthogonal with respect to another inner product.) 
Since we do not yet know what inner product to use, one might try to guess an 
inner product. Arguing by analogy, one might think that since the inner product 
in (1671) is appropriate for Hermitian Hamiltonians (if = H^), a good choice for an 
inner product associated with a PT-symmetric Hamiltonian {H = H^^) might be 

(V^,0)= f dx[i){x)f^(t){x) = [ dx[ij{-x)]*^{x), (77) 
Jc Jc 

where C is a contour in the Stokes wedges shown in Fig. [2l With this inner-product 
definition one can show by a trivial integration-by-parts argument using the time- 
independent Schrodinger equation that pairs of eigenfunctions of H associated 
with different eigenvalues are orthogonal. However, this guess for an inner product 
is not acceptable for formulating a valid quantum theory because the norm of a 
state is not necessarily positive. 

(c) The CVT inner product. To construct an inner product with a positive norm for 
a complex non-Hermitian Hamiltonian having an unbroken VT symmetry, we will 
construct a new linear operator C that commutes with both H and VT . Because C 
commutes with the Hamiltonian, it represents a symmetry oi H . We use the symbol 
C to represent this symmetry because, as we will see, the properties of C are similar 
to those of the charge conjugation operator in particle physics. The inner product 
with respect to CVT conjugation is defined as 

{m^'^^ = I rfx/^^(x)x(x), (78) 

where ip'^^'^^x) = J dyC{x, y)ip*{—y). We will show that this inner product satisfies 
the requirements for the quantum theory defined by H to have a Hilbert space with 
a positive norm and to be a unitary theory of quantum mechanics. We will represent 
the C operator as a sum over the eigenfunctions of H, but before doing so we must 
first show how to normalize these eigenfunctions. 

(d) VT -symmetric normalization of the eigenfunctions and the strange statement of 
completeness. We showed in (fT9l) that the eigenfunctions ipn{x) of H are also 
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eigenf unctions of the VT operator with eigenvalue A = e'", where A and a depend 
on n. Thus, we can construct PT-normahzed eigenf unctions 0n(x) defined by 

= e-*"/2^„(x). (79) 

By this construction, 0„(x) is still an eigenfunction of H and it is also an 
eigenfunction of VT with eigenvalue 1. One can also show both numerically and 
analytically that the algebraic sign of the VT norm in (!77|) of ipnix) is (—1)" for all 
n and for all values of e > |18] . Thus, we define the eigenfunctions so that their 
VT norms are exactly (—1)": 

/ [0„(x)]^^0„(x) = / t/x[0„(-x)]*0„(x) = (-1)", (80) 
Jc Jc 

where the contour C lies in the Stokes wedges shown in Fig. [2l In terms of these VT- 
normalized eigenfunctions there is a simple but unusual statement of completeness: 

oo 

5^(-l)>„(x)0„(y)=5(x-y). (81) 

?i=0 

This unusual statement of completeness has been verified both numerically and 
analytically to great precision for all e > [IHl EH] and a mathematical proof has 
been given [31] • It is easy to verify using flHOl) that the left side of flHTl) satisfies the 
integration rule for delta functions: / dy6{x — y)S{y — z) = 6{x — z). 

Example: VT -symmetric normalization of harmonic-oscillator eigenfunctions. For the harmonic- 
oscillator Hamiltonian H — + i^, the eigenfunctions are Gaussians multiplied by Hermite 
polynomials: ^o{x) = exp(— ^a;^), '4>i{x) = a: exp (— ^a;^) , ^'2(2;) = (2x^ — 1) exp (— ^a;^) , 
•03 (x) — [2x^ — 3a:) exp (— ix^), and so on. To normahze these eigenfunctions so that they 
are also eigenfunctions of the VT operator with eigenvalue 1, we choose 

4>q{x) = ao exp {-\x^) , 
<j)i{x) = aiia;exp {—\x'^) , 
^2{x) = a2(2a;^ - 1) exp {~\x^) , 

,/;3(a;) = a3i(2a;^ - 3a;)exp(-ia;2) , (82) 

and so on, where the real numbers a„ are chosen so that the integral in ([80)1 evaluates to (— 1)" 
for all n. It is easy to verify that if the eigenfunctions (j)n{x) are substituted into ([81]) and the 
summation is performed, then the result is the Dirac delta function on the right side of ((8T|l . 

(e) Coordinate-space representation of H and G and the spectral Zeta function. 
From the statement of completeness in ( IHTi) we can construct coordinate-space 
representations of the linear operators. For example, since the coordinate-space 

representation of the parity operator is V{x,y) = 6{x + y), we have 

00 

v{x,y) = J2{-iy^Mx)M-y)- (83) 

n=0 

We can also construct the coordinate-space representations of the Hamiltonian and 
the Green's function, 

00 

H{x,y)= J2i-'^TEnMx)My), 

n=0 
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oo 

G{x,y) = ^(-1)"— (84) 

n=0 

and using (IHOj) it is straightforward to show that G is the functional inverse of H: 
J dy H{x,y)G{y, z) = 6{x — z). For the class of PT-symmetric Hamiltonians in 
(|T2|) this equation takes the form of a differential equation satisfied by G(x, y): 

+ x^itxY^ Gix, y) = 6ix - y). (85) 

Equation ( !85l) can be solved in terms of associated Bessel functions in each of two 
regions, x > y and x < y. The solutions can then be patched together ed, x = y to 
obtain a closed-form expression for G{x,y) [HJ [50]. One then uses (!75|) to find an 
exact formula for the spectral zeta function for all values of e: 



^ E 



\ ^ cos (^) sin (^) 



(4 + e)^r(i±f)r(|±f) 



(86) 



(f) Construction of the C operator. The situation here in which half of the energy 
eigenstates have positive norm and half have negative norm [see ( IHOj) ] is analogous 
to the problem that Dirac encountered in formulating the spinor wave equation 
in relativistic quantum theory ^1]. Following Dirac, we attack the problem of 
an indefinite norm by finding an interpretation of the negative-norm states. For 
any Hamiltonian H having an unbroken VT symmetry there exists an additional 
symmetry of H connected with the fact that there are equal numbers of positive- 
and negative-norm states. The linear operator C that embodies this symmetry can 
be represented in coordinate space as a sum over the PT-normalized eigenfunctions 
of the PT-symmetric Hamiltonian in ( |T2l) : 

oo 

c{x,y) = J2M^)My)- (87) 

n=0 

Notice that this equation is identical to the statement of completeness in f lHTj) except 
that the factor of (—1)" is absent. We can use ( IHOi) and ( IHTi) to verify that the square 
of C is unity (C^ = 1): 

dyC{x,y)C{y,z) = 6{x - z). (88) 

Thus, the eigenvalues of C are ±1. Also, C commutes with H. Therefore, since C 
is linear, the eigenstates of H have definite values of C. Specifically, 



C(j)n{x) = J dyC{x,y)(j)n{y) 

oo „ 

= ^0™(a;) / dy<j)m{y)My) = {-ITM^)- (89) 

This new operator C resembles the charge-conjugation operator in quantum field 
theory. However, the precise meaning of C is that it represents the measurement of 
the sign of the VT norm in (l80ll of an eigenstate. 
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The operators V and C are distinct square roots of the unity operator 6{x — y). 
That is, = = 1, but V ^ C because V is real, while C is complex. The 
parity operator in coordinate space is explicitly real (V{x,y) = 6{x + y)), while the 
operator C{x,y) is complex because it is a sum of products of complex functions. 
The two operators V and C do not commute. However, C does commute with VT. 

(g) Positive norm and unitarity in VT -symmetric quantum mechanics. Having 
constructed the operator C, we can now use the new CVT inner-product defined 
in (1781) . Like the VT inner product, this new inner product is phase independent. 
Also, because the time-evolution operator (as in ordinary quantum mechanics) is 
e~*^* and because H commutes with VT and with CVT, both the VT inner product 
and the CVT inner product remain time independent as the states evolve. However, 
unlike the VT inner product, the CVT inner product is positive definite because C 
contributes a factor of —1 when it acts on states with negative VT norm. In terms 
of the CVT conjugate, the completeness condition reads 



4-3. Comparison of Hermitian and VT- Symmetric Quantum Theories 

We have shown in Subsecs. 14.11 and 14.21 how to construct quantum theories based on 
Hermitian and on non-Hermitian Hamiltonians. In the formulation of a conventional 
quantum theory defined by a Hermitian Hamiltonian, the Hilbert space of physical states 
is specified even before the Hamiltonian is known. The inner product ( 1671) in this vector 
space is defined with respect to Dirac Hermitian conjugation (complex conjugate and 
transpose). The Hamiltonian is then chosen and the eigenvectors and eigenvalues of 
the Hamiltonian are determined. In contrast, the inner product for a quantum theory 
defined by a non-Hermitian PT-symmetric Hamiltonian depends on the Hamiltonian 
itself and thus is determined dynamically. One must solve for the eigenstates of H before 
knowing the Hilbert space and the associated inner product of the theory because the 
C operator is defined and constructed in terms of the eigenstates of the Hamiltonian. 
The Hilbert space, which consists of all complex linear combinations of the eigenstates 
of H, and the CVT inner product are determined by these eigenstates. 

The operator C does not exist as a distinct entity in ordinary Hermitian quantum 
mechanics. Indeed, if we allow the parameter e in (1121) to tend to 0, the operator C in this 
limit becomes identical to V and the CVT operator becomes T, which performs complex 
conjugation. Hence, the inner product defined with respect to CVT conjugation reduces 
to the inner product of conventional quantum mechanics and ( ISTI) reduces to the usual 
statement of completeness Xln '^"(^)'^n(?/) = ^(^ ~ v)- 

The CVT inner product is independent of the choice of integration contour C as 
long as C lies inside the asymptotic wedges associated with the boundary conditions 
for the eigenvalue problem. In ordinary quantum mechanics, where the positive-definite 
inner product has the form J dx f*{x)g{x), the integral must be taken along the real 



oo 




(90) 



n=0 
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axis and the path of integration cannot be deformed into the complex plane because the 
integrand is not analytic. The VT inner product shares with the CVT inner-product 
the advantage of analyticity and path independence, but it suffers from nonpositivity. It 
is surprising that we can construct a positive-definite metric by using CVT conjugation 
without disturbing the path independence of the inner-product integral. 

Time evolution is expressed by the operator e~*^* whether the theory is determined 
by a PT-symmetric Hamiltonian or just an ordinary Hermitian Hamiltonian. To 
establish unitarity we must show that as a state vector evolves, its norm does not change 
in time. If iPq{x) is any given initial vector belonging to the Hilbert space spanned 
by the energy eigenstates, then it evolves into the state iptix) at time t according to 
ipt{x) = e'^^^ipQ^x). With respect to the CVT inner product the norm of ipt{x) does 
not change in time because H commutes with CVT . 

4.4- Observables 

How do we represent an observable in "PT-symmetric quantum mechanics? Recall 
that in ordinary quantum mechanics the condition for a linear operator A to be an 
observable is that A = This condition guarantees that the expectation value of 
A in a state is real. Operators in the Heisenberg picture evolve in time according to 
A{t) = e*^*A(0)e~*^*, so this Hermiticity condition is maintained in time. In PT- 
symmetric quantum mechanics the equivalent condition is that at time t = the 
operator A must obey the condition A^ = CVTACVT, where A^ is the transpose of A 
|48] . If this condition holds at t = 0, then it will continue to hold for all time because 
we have assumed that H is symmetric {H = H^). This condition also guarantees that 
the expectation value of A in any state is realj^ 

The operator C itself satisfies this requirement, so it is an observable. The 
Hamiltonian is also an observable. However, the x and p operators are not observables. 
Indeed, the expectation value of x in the ground state is a negative imaginary number. 
Thus, there is no position operator in PT-symmetric quantum mechanics. In this sense 
PT-symmetric quantum mechanics is similar to fermionic quantum field theories. In 
such theories the fermion field corresponds to the x operator. The fermion field is 
complex and does not have a classical limit. One cannot measure the position of an 
electron; one can only measure the position of the charge or the energy of the electron! 

One can see why the expectation of the x operator in PT-symmetric quantum 
mechanics is a negative imaginary number by examining a classical trajectory like that 
shown in Fig. [151 Note that this classical trajectory has left-right (PT) symmetry, 
but not up-down symmetry. Also, the classical paths favor (spend more time in) the 
lower-half complex-x plane. Thus, the average classical position is a negative imaginary 
number. Just as the classical particle moves about in the complex plane, the quantum 

tt The requirement given here for A to be an observable involves matrix transposition. This condition 
is more restrictive than is necessary and it has been generalized by Mostafazadeh. See Refs. [52 l [53 l [54] . 
Note that if the matrix transpose symmetry condition on the Hamiltonian is removed we must introduce 
a biorthogonal basis. See Refs. [31] and [47] . 
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probability current flows about in the complex plane. It may be that the correct 
interpretation is to view PT-symmetric quantum mechanics as describing the interaction 
of extended, rather than pointlike objects. 

4-5. Pseudo-Hermiticity andVT Symmetry 

The thesis of this paper - replacing the mathematical condition of Hermiticity by the 
more physical condition of VT symmetry - can be placed in a more general mathematical 
context known as pseudo-Hermiticity. A linear operator A is pseudo-Hermitian if there 
is a Hermitian operator rj such that 



The operator rj is often called an intertwining operator. The condition in fl9T|) obviously 
reduces to ordinary Hermiticity when the intertwining operator t] is the identity 1. The 
concept of pseudo-Hermiticity was introduced in the 1940s by Dirac and Pauli, and later 
discussed by Lee, Wick, and Sudarshan, who were trying to resolve the problems that 
arise in quantizing electrodynamics and other quantum field theories in which negative- 
norm states states appear as a consequence of renormalization [551 EHl EZl EH [59] . These 
problems are illustrated very clearly by the Lee model, which is discussed in Subsec. 18.31 
Mostafazadeh observed that because the parity operator V is Hermitian, it may 
be used as an intertwining operator. The class of Hamiltonians H in ( fT2l) is pseudo- 
Hermitian because the parity operator V changes the sign of x while Dirac Hermitian 
conjugation changes the sign of i [60l [611 [62l [53] : 



For some other references on the generalization of VT symmetry to pseudo-Hermiticity 
see [63] . 

5. Illustrative 2x2 Matrix Example of a PT-Symmetric Hamiltonian 

It is always useful to study exactly solvable models when one is trying to understand 
a formal procedure like that discussed in Sec. [H One exactly solvable model, which 
has been used in many papers on VT symmetry, is due to Swanson [64]. This model is 
exactly solvable because it is quadratic in x and p. In this section we use an even simpler 
model to illustrate the construction of a quantum theory described by a PT-symmetric 
Hamiltonian. We consider the elementary 2x2 Hamiltonian matrix 



where the three parameters r, s, and 6 are real [65J. The Hamiltonian in (l93ll is not 
Hermitian, but it is easy to see that it is VT symmetric, where we define the parity 
operator as 



= r]Ar]. 



(91) 



= VHV. 



(92) 




(93) 




(94) 
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and we define the operator T to perform complex conjugation. 

As a first step in analyzing the Hamiltonian fl93l) . we calculate its two eigenvalues: 

E± = rcose±{s^ -r'^sin^ef'^. (95) 

There are clearly two parametric regions to consider, one for which the square root in 
(!95|) is real and the other for which it is imaginary. When < sin^ 9, the energy 
eigenvalues form a complex conjugate pair. This is the region of broken VT symmetry. 
On the other hand, when > sin^ 6', the eigenvalues e±= r cos ^ ± (s^ — sin^ 6)^^"^ 
are real. This is the region of unbroken VT symmetry. In the unbroken region the 
simultaneous eigenstates of the operators H and VT are 

1 / pW2 \ i ( 

= I and \E.) = I , (96) 



cos« 



where 



r 



sin a = - sin 6. (97) 

s 

The VT inner product gives 

(E±,E±)=±1 and (E±,^^) = 0, (98) 

where (m, v) = iVTu) ■ v. With respect to the VT inner product, the vector space 
spanned by the energy eigenstates has a metric of signature (+,—). If the condition 
> sin^ 9 for an unbroken VT symmetry is violated, the states ( !96l) are no longer 
eigenstates of VT because a becomes imaginary. When VT symmetry is broken, the 
VT norm of the energy eigenstate vanishes. 

Next, we construct the operator C using fl57|) : 




C = , . . . (99) 

cos a 

Note that C is distinct from H and V and it has the key property that 

C\E±) = ±\E±). (100) 

The operator C commutes with H and satisfies = 1. The eigenvalues of C are precisely 
the signs of the VT norms of the corresponding eigenstates. Using the operator C we 
construct the new inner product structure {u\v) = (CVTu) ■ v. This inner product 
is positive definite because {E±\E±) = 1. Thus, the two-dimensional Hilbert space 
spanned by \E±), with inner product (-I-), has signature (+,+). 

Finally, we prove that the CVT norm of any vector is positive. For the arbitrary 
vector ip = (^), where a and b are any complex numbers, we see that 

= f f V -PTi: 

cvT4,=^ I ^;+^^;^^"" ). (101) 

cos a \ b — la sm a 
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Thus, = {CVT ip)-ip = -^[a*a + b*b + i{b*b- a*a) sin a]. Now let a = x + iy and 

b = u + iv, where x, y, u, and v are real. Then 

(i'li') = (x^ + f ^ + 2xvsma + y'^ + — 2yusma) , (102) 

cos a 

which is exphcitly positive and vanishes only ifx = y = u = v = 0. 

Since {u\ denotes the CPT-conjugate of \u), the completeness condition reads 



E+){E+\ + \E^){E^ 




(103) 



Furthermore, using the CVT conjugate {E±\, we can represent C as 

C = \E+){E+\-\E^){E^\. (104) 

In the limit 6^0, the Hamiltonian (193!) for this two-state system becomes 
Hermitian and C reduces to the parity operator V. Thus, CVT invariance reduces 
to the standard condition of Hermiticity for a symmetric matrix; namely, H = H*. 



6. Calculation of the C Operator in Quantum Mechanics 

The distinguishing feature of "PT-symmetric quantum mechanics is the C operator. In 
ordinary Hermitian quantum mechanics there is no such operator. Only a non- Hermitian 
PT-symmetric Hamiltonian possesses a C operator distinct from the parity operator V. 
Indeed, if we were to sum the series in fl87|) for a "PT-symmetric Hermitian Hamiltonian, 
the result would be P, which in coordinate space is 6{x + y). [See (l83l).] 

While the C operator is represented formally in (IHTl) as an infinite series, it is not 
easy to evaluate the sum of this series. Calculating C by direct brute-force evaluation 
of the sum in fl87l) is not easy in quantum mechanics because it is necessary to find 
all the eigenfunctions (pnix) of H. Furthermore, such a procedure cannot be used in 
quantum field theory because in field theory there is no simple analog of the Schrodinger 
eigenvalue differential equation and its associated coordinate-space eigenfunctions. 

The first attempt to calculate C relied on a perturbative approach [66] ■ In this 
paper the PT-symmetric Hamiltonian 

H = + \x^ + itx^ (105) 

is considered, where e is treated as a small real parameter. When e = 0, the 
Hamiltonian reduces to the Hamiltonian for the quantum harmonic oscillator, all of 
whose eigenfunctions can be calculated exactly. Thus, it is possible to express the 
eigenfunctions of H in (11051) in the form of perturbation series in powers of e. For each 
of the eigenfunctions the first few terms in the perturbation series were calculated. 
These perturbation series were then substituted into flHTl) and the summation over 
the nth eigenfunction was performed. The result is a perturbation-series expansion 
of the C operator. This calculation is long and tedious and the final result is quite 
complicated. However, the calculation turned out to be of great value because while 
the final answer is complicated, it was discovered that the answer in coordinate space 
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simplified dramatically if the C operator is written as the exponential of a derivative 
operator Q multiplying the parity operator V: 



C{x,y) = exp 



6{x + y). (106) 



The simplification that occurs when C is written in the form C = e^V is that while the 
expression for C is a series in all positive integer powers of e, Q is a. series in odd powers 
of e only. Since Q is a series in odd powers of e, in the limit e — > the function Q 
vanishes. Thus, in this limit the C operator tends to the parity operator V. 

The expression in fll06p need not be limited to coordinate space. A more general 
way to represent the C operator is to express it generically in terms of the fundamental 
dynamical operators ct and p: 

C = e^(^'P)p. (107) 

Written in this form, Q is a real function of its two variables. By seeking the C operator 
in the form (11071) we will be able to devise powerful analytic tools for calculating it. 

We illustrate the representation C = e^V by using two elementary Hamiltonians. 
First, consider the shifted harmonic oscillator H = + + iex. This Hamiltonian 
has an unbroken VT symmetry for all real e. Its eigenvalues En = n + ^ + ^e^ are all 
real. The exact formula for C for this theory is given exactly by C = e^V, where 

Q = -ep. (108) 

Note that in the limit e — > 0, where the Hamiltonian becomes Hermitian, C becomes 
identical with V. 

As a second example, consider the non-Hermitian 2x2 matrix Hamiltonian (J93l) . 
The associated C operator in fl99l) can be rewritten in the form C = e'^V, where 

e = (109) 

^ \ 1 + sm a y 

and 

^2 = I : "1 (110) 




is the Pauli sigma matrix. Again, observe that in the limit 6^0, where the Hamiltonian 
becomes Hermitian, the C operator becomes identical with V. 



6.1. Algebraic Equations Satisfied by the C Operator 

Fortunately, there is a relatively easy algebraic way to calculate the C operator, and the 
procedure circumvents the difficult problem of evaluating the sum in ( 1571) . As a result, 
the technique readily generalizes from quantum mechanics to quantum field theory. In 
this subsection we show how to use this technique to calculate C for the PT-symmetric 
Hamiltonian in (11051) . We explain how to calculate C perturbatively to high order in 
powers of e for this cubic Hamiltonian. Calculating C for other kinds of interactions is 
more difficult and may require the use of semiclassical approximations [67] . 
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To calculate C we make use of its three crucial algebraic properties. First, C 
commutes with the space-time reflection operator VT, 

[c,vr] = o, (111) 

although C does not, in general, commute with V or T separately. Second, the square 
of C is the identity, 

C^ = l, (112) 

which allows us to interpret C as a reflection operator. Third, C commutes with H, 

[C,H]=0, (113) 

and thus is time independent. To summarize, C is a time-independent PT-symmetric 
reflection operator. 

The procedure for calculating C is simply to substitute the generic operator 
representation in (11071) into the three algebraic equations (lllip - (I113P in turn and 
to solve the resulting equations for the function Q. First, we substitute (I107p into the 
condition (II lip to obtain 

^Q(x,p) _ 'pq-^Q{x,p)'pq- _ ^Q(-x,p) ^ (114) 

from which we conclude that Q{x,p) is an even function of x. 

Second, we substitute (I107P into the condition (11120 and find that 

^Qix,p)'p^Q{x,p)jy _ ^Q{x,p) ^Q{-x~p) _ (115) 

which implies that Q{x,p) = —Q{—x, —p). Since we already know that Q{x,p) is an 
even function of x, we conclude that it is also an odd function of p. 

The remaining condition (IllSp to be imposed is that the operator C commutes 
with H. While the first two conditions are, in effect, kinematic conditions on Q that are 
generally true for any Hamiltonian, condition (I113p is equivalent to imposing the specific 
dynamics of the particular Hamiltonian that defines the quantum theory. Substituting 
C = e^^'^'P^V into ffTT3D . we get 

eQ(^-p)[p,if] + [e«(^'P),ff]p = 0. (116) 

This equation is difficult to solve in general, and to do so we must use perturbative 
methods, as we explain in the next subsection. 

6.2. Perturbative Calculation ofC. 

To solve (I116p for the Hamiltonian in (llOSp . we express this Hamiltonian in the form 
H = Hq + eHi. Here, Hq is the harmonic-oscillator Hamiltonian Hq = \fP' + ^x^, which 
commutes with the parity operator P, and Hi = ix^, which anticommutes with V. 
Thus, the condition (I116P becomes 

2ge«(^'P)i7i = [e«(^'P),/7]. (117) 

We expand the operator Q{x,p) as a perturbation series in odd powers of e: 

Q{x,p) = eQi{x,p) + e^Q3{x,p) + e'Q^ix^p) + ■ ■ ■ . (118) 
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Substituting the expansion in flllSp into the exponential e*5(^'P), we get after some 
algebra a sequence of equations that can be solved systematically for the operator- 
valued functions Qn{x,p) (n = 1,3,5,...) subject to the symmetry constraints that 
ensure the conditions fillip and flll2p . The first three of these equations are 

[Hq, Qi] = — 2Hi, 

[Ho,Q3] = -l[Qi,[Qi,H,]], 

[Ho,Q5] = [^1' [^1' [Qi' lQl^H^]]]] - |[Qi, [Qs,H^]] 

+ l[Q,,[Q,,H,]]. (119) 

Let us solve these equations for the Hamiltonian in (11051) . for which Hq = |]3^ + 
and Hi = ix^ . The procedure is to substitute the most general polynomial form for 
Qn using arbitrary coefficients and then to solve for these coefficients. For example, to 
solve [Hq, Qi] = —2ix^, the first of the equations in (11191) . we take as an ansatz for Qi 
the most general Hermitian cubic polynomial that is even in x and odd in p\ 

Qiix,p) = Mf + Nxpx, (120) 

where M and are numerical coefficients to be determined. The operator equation for 
Qi is satisfied if M = — | and = —2. 

It is algebraically tedious but completely straightforward to continue this process. 
In order to present the solutions for Qn{x,p) {n > 1), it is convenient to introduce 
the following notation: Let Sm,n represent the totally symmetrized sum over all terms 
containing m factors of p and n factors of x. For example, 

5*0,0 = 1, 5*0,3 = x^, 5*1,1 = ^ixp + px) , 5*1,2 = \ [x^P + xpx + px^) , (121) 

and so on. (The properties of the operators Sm,n are summarized in Ref. |68j.) 

In terms of the symmetrized operators S^.n the first three functions Q2n+i are 

Qi = - - 2^1,2, 

^3 = Wp' + f^3,2 + 85i,4-12p, 

Q5 = - W - ^^5,2 - ^^3,4 - 645i,6 + + ^^1,2. (122) 

Together, (I107p . (IllSp . and (11221) represent an explicit perturbative expansion of C in 
terms of the operators x and p, correct to order e^. 

To summarize, using the ansatz (llOTp we can calculate C to high order in 
perturbation theory. We are able to perform this calculation because this ansatz obviates 
the necessity of calculating the PT-normalized wave functions 0„(a;). We show how use 
these same techniques for quantum field theory in Sec. [HI 

6. 3. Perturbative Calculation of C for Other Quantum- Mechanical Hamiltonians 

The C operator has been calculated perturbatively for a variety of quantum-mechanical 
models. For example, let us consider first the case of the Hamiltonian 

H = \{p' + q^) +\ {x^ + y^) + lex^y, (123) 
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which has two degrees of freedom. The energy levels of this complex Henon-Heiles 
theory were studied in Ref. p9]. The C operator for this Hamiltonian was calculated in 
[Tot EH E2]. The perturbative result for the Q = Qie + Q^e^ + . . . operator is 

Qi (x, y, PA)= - ^p^Q. - \Si^iy - I^^q, 

Q3ix,y,p,q) = IIpY + Ip^g + ^5i,iT2,i 

256 -2^ I 512 Q I 288 r. ^ 32 - 2^53 



405^ J^'^ ' 405 •^'^y ' 405^'^^ 405 
,736 -2^ 256 r. ^-.S 1 608 o s, 128-4- 8^ (^0A\ 

where T^^n represents a totally symmetric product of m factors of q and n factors of y. 
For the Hamiltonian 

H = \{p^ + +f'^) +\ [x^ + f + z^) + iexyz, (125) 
which has three degrees of freedom, we have [70l [71], [72] 

Qi{x,y, z,p,q,f) = - ^{yzp + xzq + xyf) - ^pqf, 

Q3{x, y, z, p, g, f) = III [p^qr + fpf + r^qp) 

+ ^\pxp{yf + zq) + qyq{xf + zp) + rzr{xq + yp)] 
— {xpxqr + yqypr + zfzpq) + |^ (xpxyz + yqyxz 
+ zfzxy) — -^[x^iyr + zq) + y^{xr + zp) 
+ z^{xq + yp)] - {p^yz + q^xz + r^xy) . (126) 

An extremely interesting and deceptively simple quantum-mechanical model is the 
PT-symmetric square well, whose Hamiltonian on the domain < a; < vr is given by 

H = p^ + V{x), (127) 

where in the coordinate representation V{x) = oo for x < and x > n and 

^ ^ ^ -ze for < X < f . ^ ' 



The PT-symmetric square-well Hamiltonian was invented by Znojil [73] and it has been 
heavily studied by many researchers [TH, [TSj, [S21 [TH] • This Hamiltonian reduces to the 
conventional Hermitian square well in the limit as e — 0. For H in fll27p the parity 
operator V performs a reflection about x = ^: V : x ^ — x. 

In all the examples discussed so far the coordinate-space representation of the C 
operator is a combination of integer powers of x and integer numbers of derivatives 
multiplying the parity operator V. Hence, the Q operator is a polynomial in the 
operators x and p = —i^- The novelty of the PT-symmetric square- well Hamiltonian 
(11271) and (11281) is that C contains integrals of V. Thus, the Q operator, while it is a 
simple function, is not a polynomial in x and p and therefore it cannot be found easily 
by the algebraic perturbative methods introduced in Subsec. 16.11 

Instead, in Ref. [77] the C operator for the square-well Hamiltonian was calculated 
by using the brute-force approach of calculating the PT-normalized eigenfunctions 
0n(x) of H and summing over these eigenfunctions. The eigenfunctions (()n{x) were 



Making Sense of Non-Hermitian Hamiltonians 



49 



obtained as perturbation series to second order in powers of e. The eigenfunctions were 
then normahzed according to (|80|) . substituted into the formula fl87|) . and the sum was 
evaluated directly to obtain the C operator accurate to order e^. The advantage of the 
domain < x < vr is that this sum reduces to a set of Fourier sine and cosine series that 
can be evaluated in closed form. After evaluating the sum, the result was translated to 
the symmetric region — | < a; < ^. On this domain the parity operator in coordinate 
space is V{x, y) = 6{x + y). 

The last step in the calculation is to show that the C operator to order has the 
form e'^V and then to evaluate the function Q to order e^. The final result for Q{x,y) 
on the domain — | < a; < ^ has the relatively simple structure 

Q{x, y) = \it[x -y + e{x-y){\x + y\-Ti)] + O^e^), (129) 

where e{x) is the standard step function 

r 1 (x>o), 

e{x) = I (x = 0), (130) 
[ -1 (x < 0). 

When we say that the formula for Q{x,y) has a relatively simple structure, we 
mean that this structure is simple in comparison with the formula for the C operator 
expressed as an expansion C{x, y) = C^^^ + eC^^^ + e^C^^^ + 0{e^). The formulas for the 
first three coefficients in this series for C are 

C^^\x,y) = 6{x + y), 

C^^\x,y) = \i[x + y + e{x + y) i\x - y\ -vr)], 

C(2)(x, y) = ^vr^ - ^(x^ + y^) ^(x + y) - j-^{y^ - x^) e{y - x) 

+ |x|/7r - j^7T^{x + y) e{x + y) + |7r(x|x| + y\y\) e{x + y) 

- \xy{\x\[e{x - y) e{-x - y) + e{y - x) e{x + y)] 

+ \y\\e{y - X) ei-x -y) + Oix - y) e{x + y)]}. (131) 

We plot the imaginary part of C^^-'(x,?/) in Fig. [22], and we plot C^'^\x,y) in Fig. [231 
These three-dimensional plots show C^^\x,y) and C^'^\x,y) on the symmetric domain 

-f < {x,y) < f. 

The most noteworthy property of the C operator for the square-well model is that 
the associated operator Q is a nonpolynomial function, and this kind of structure had 
not been seen in previous studies of C. It was originally believed that for such a simple 
PT-symmetric Hamiltonian it would be possible to calculate the C operator in closed 
form. It is a surprise that for this elementary model the C operator is so nontrivial. 

6.4- Mapping from a VT -Symmetric Hamiltonian to a Hermitian Hamiltonian 

Mostafazadeh observed that the square root of the positive operator can be used 
to construct a similarity transformation that maps a non-Hermitian PT-symmetric 
Hamiltonian H to an equivalent Hermitian Hamiltonian h [60l [78] : 

h = e-«/2^e'3/2. (132) 
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Figure 22. Three-dimensional plot of the imaginary part of C^^\x,y), the first-order 
pcrtm'bative contribution in (|13ip to the C operator in coordinate space. The plot is 
on the symmetric square domain < {x,y) < ^. Note that C^^'>{x,y) vanishes on 
the boundary of this square domain because the eigenfunctions (j)n{x) are required to 
vanish at a; = and x = n. 




Figure 23. Three-dimensional plot of C^'^\x,y) in ()13ip on the symmetric square 
domain — < {x,y) < ^. The function C'-^^(a;,y) vanishes on the boundary of this 
square domain because the eigenfunctions (j)n (x) from which it was constructed vanish 
at the boundaries of the square well. 

He noted that h is equivalent to H because it has the same eigenvalues as H. 

To understand why (11321) is valid, recall from f llOTp that the C operator has the 
general form C = e^V, where Q = Q{x,p) is a Hermitian function of the fundamental 
dynamical operator variables of the quantum theory. Multiplying C on the right by V 
gives an expression for e*^: 

= CV, (133) 

which indicates that CV is a positive and invertible operator. 
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To verify that h in (11321) is Hermitian, take the Hermitian conjugate of (11321) 

= e«/2i7te-Q/2^ (134) 

which can be rewritten as 

Use ffT33D to replace by CV and e'^ by VC 

h) = e-^/^CVH^VCe^/^ (136) 
and recall from (l92l) that H"^ can be replaced by VHV. This gives 

/it = e-^'^CVVHVVCe^'^ = e-^'^CHCe^'^ . (137) 

Finally, one use the fact that C commutes with H [see ( 1113P ] and that the square of C is 
unity [see (11121) ] to reduce the right side of (11361) to the right side of (11321) . This verifies 
that h is Hermitian in the Dirac sense. 

We conclude from this calculation that for every non-Hermitian "PT-symmetric 
Hamiltonian H whose VT symmetry is unbroken, it is, in principle, possible to construct 
via (11321) a Hermitian Hamiltonian h that has exactly the same eigenvalues as H. Note 
that it is crucial that H have an unbroken VT symmetry because this allows us to 
construct C, which in turn allows us to construct the similarity operator e*^/^. 

This construction poses the following question: Are PT-symmetric Hamiltonians 
physically new and distinct from ordinary Hermitian Hamiltonians, or do they describe 
exactly the same physical processes that ordinary Hermitian Hamiltonians describe? 

There are two answers to this question, the first being technical and practical and 
the second being an answer in principle. First, while it is theoretically possible to 
construct the Hermitian Hamiltonian h whose spectrum is identical to that of if, it 
cannot in general be done except at the perturbative level because the transformation 
is so horribly complicated. (See, for example, the discussion of the square well in 
Subsec. 16.31 ) Furthermore, while the PT-symmetric Hamiltonian H is simple in 
structure and easy to use in calculations because the interaction term is local, it is 
shown in Ref. [79] that h is nonlocal (its interaction term has arbitrarily high powers 
of the variables x and p.) Thus, it is not just difficult to calculate using h] it is almost 
impossible because the usual regulation schemes are hopelessly inadequate. 

There is only one known nontrivial example for which there is actually a closed-form 
expression for both H and h, and this is the case of the quartic Hamiltonian discussed 
in Subecs. 12.61 - 12.81 Even this case, it is not possible to construct the C operator in 
closed form because this operator is nonlocal (it contains a Fourier transform) and it 
performs a transformation in the complex plane. This is the only known example for 
which it is practical to calculate with both H and h. Hence, while the mapping from H 
to h is of great theoretical interest, it does not have much practical value. 

The second answer is of greater importance because it leads immediately to physical 
considerations. The transformation from if to /i in (11321) is a similarity and not a unitary 
transformation. Thus, while the eigenvalues of H and h are the same, relationships 
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between vectors are changed; pairs of vectors that are orthogonal are mapped into 
pairs of vectors that are not orthogonal. Thus, experiments that measure vectorial 
relationships can distinguish between H and h. One plausible experiment, which 
involves the speed of unitary time evolution, is described in detail in Subsec. 17.31 

7. Applications of PT-Symmetric Hamiltonians in Quantum Mechanics 

It is not yet known whether non-Hermitian PT-symmetric Hamiltonians describe 
phenomena that can be observed experimentally. However, non-Hermitian VT- 
symmetric Hamiltonians have already appeared in the literature very often and their 
remarkable properties have been noticed and used by many authors. For example, in 
1959 Wu showed that the ground state of a Bose system of hard spheres is described by a 
non-Hermitian Hamiltonian |80]. Wu found that the ground-state energy of this system 
is real and he conjectured that all the energy levels were real. HoUowood showed that the 
non-Hermitian Hamiltonian of a complex Toda lattice has real energy levels [8]. Cubic 
non-Hermitian Hamiltonians of the form H = + ix^ (and also cubic quantum field 
theories) arise in studies of the Lee- Yang edge singularity [HI |82[ [831 El] and in various 
Reggeon field theory models [85l[86]. In all of these cases a non-Hermitian Hamiltonian 
having a real spectrum appeared mysterious at the time, but now the explanation is 
simple: In every case the non-Hermitian Hamiltonian is VT symmetric. In each case 
the Hamiltonian is constructed so that the position operator x or the field operator 
(f) is always multiplied by i. Hamiltonians having VT symmetry have also been used 
to describe magnetohydrodynamic systems [HTJ EH] and to study nondissipative time- 
dependent systems interacting with electromagnetic fields {89] . 




In this section we describe briefly four areas of quantum mechanics in which non- 
Hermitian "PT-Hamiltonians play a crucial and significant role. 

7.1. New VT- Symmetric Quasi-Exactly Solvable Hamiltonians 

A quantum-mechanical Hamiltonian is said to be quasi-exactly solvable (QES) if a finite 
portion of its energy spectrum and associated eigenfunctions can be found exactly and 
in closed form ^0]. QES potentials depend on an integer parameter J; for positive 
values of J one can find exactly the first J eigenvalues and eigenfunctions, typically of 
a given parity. QES systems are classified using an algebraic approach in which the 
Hamiltonian is expressed in terms of the generators of a Lie algebra ^92j . This approach 
generalizes the dynamical-symmetry analysis of exactly solvable quantum-mechanical 
systems whose entire spectrum may be found in closed form by algebraic means. 

Prior to the discovery of non-Hermitian PT-symmetric Hamiltonians the lowest- 
degree one- dimensional QES polynomial potential that was known was a sextic potential 
having one continuous parameter as well as a discrete parameter J. A simple case of 
such a potential is [91] 




V{x) = X 



.6 



(4J- l)xl 



(138) 
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For this potential, the Schrodinger equation —ip"{x) + \V{x) — E\'4){x) = has J even- 
parity solutions of the form 

j-i 

^(x) = e-"*/^ ^ CfcX^^ (139) 

fc=0 

The coefficients for < A; < J — 1 satisfy the recursion relation 

4(J - k)ck-i + Eck + 2{k + l){2k + l)ck+i = 0, (140) 

where c_i = cj = 0. The linear equations (I140p have a nontrivial solution for 
Co, Ci, cj-i if the determinant of the coefficients vanishes. For each integer J this 
determinant is a polynomial of degree J in the variable E. The roots of this polynomial 
are all real and they are the J quasi-exact energy eigenvalues of the potential fll38p . 

The discovery of VT symmetry allows us to introduce an entirely new class of 
QES quartic polynomial potentials having two continuous parameters in addition to the 
discrete parameter J. The Hamiltonian has the form [93] 

H = p^ -x^ + 2iax^ + (a^ - 2b)x^ + 2i{ab - J)x, (141) 

where a and b are real and J is a positive integer. The spectra of this family of 
Hamiltonians are real, discrete, and bounded below. Like the eigenvalues of the potential 
(11381) . the lowest J eigenvalues of H are the roots of a polynomial of degree J. 

The eigenfunction tlj{x) satisfies "PT-symmetric boundary conditions; it vanishes in 
the Stokes wedges shown in Fig. |H The eigenfunction satisfies 

-^"(x) + [-x^ + 2iax^ + (a^ - 2b)x^ + 2i{ab - J)x] ^{x) = Eip{x).{U2) 

We obtain the QES portion of the spectrum of H in ( 114ip by making the ansatz 
ip{x) = exp (— — ^ax^ — ibx) Pj-i{x), where 

J-2 

Pj.i{x)=x''' + J2ckx'' (143) 

A:=0 

is a polynomial in x of degree J — 1. Substituting ip{x) into the differential equation 
(11421) . dividing off the exponential, and collecting powers of x, we obtain a polynomial 
in X of degree J — 1. Setting the coefficients of x'^ (1 < < J — 1) to gives a system 
of J — 1 simultaneous linear equations for the coefficients {0 < k < J — 2). We solve 
these equations and substitute the values of Ck into the coefficient of x°. This gives a 
polynomial Qj{E) of degree J in the eigenvalue E. The coefficients of this polynomial 
are functions of the parameters a and b. The first two of these polynomials are 

Q-^ = E - b'^ - a, 

Q^ = - (26^ + 4a)E + b^ + Aab'^ - 46 + 30^. (144) 

The roots of Qj{E) are the QES portion of the spectrum of H. 

The polynomials Qj{E) simplify dramatically when we substitute 

E = E + b^ + Ja and K = 4b + a^. (145) 
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The new polynomials then have the form 



Qi 
Q2 
Q3 
Q4 
Q5 



F^ 
F' 



AKF - 16, 

lOKF'^ - 96F + 9K'^, 

20KF^ - 336^2 + GiK^F + 768K. 



K. 



(146) 



The roots of these polynomials are all real as long as > -ft'criticah where -ft'criticai is a 
function of J. 

Had PT-symmetric quantum mechanics not been discovered, this beautiful family 
of quartic QES Hamiltonians would never have been considered because the quartic 
term has a negative sign. Lacking the discovery that PT-symmetric Hamiltonians have 
positive spectra, the Hamiltonians in (11411) would have been rejected because it would 
have been assumed that the spectra of such Hamiltonians would be unbounded below. 
Quasi-exactly solvable sextic PT-symmetric Hamiltonians are studied in [9^. 

7.2. Complex Crystals 

An experimental signal of a complex Hamiltonian might be found in the context of 
condensed matter physics. Consider the complex crystal lattice whose potential is 
V{x) = ism X. The optical properties of complex crystal lattices were first studied 
by Berry and O'Dell, who referred to them as complex diffraction gratings [95] . 




While the Hamiltonian H = + i sin x is not Hermitian, it is VT symmetric and 
all of its energy bands are real. At the edge of the bands the wave function of a particle 
in such a lattice is bosonic (27r-periodic), and unlike the case of ordinary crystal lattices, 
the wave function is never fermionic (47r-periodic). The discriminant for a Hermitian 
sin(x) potential is plotted in Fig. [2l]and the discriminant for a non-Hermitian ism{x) 
potential is plotted in Fig. [25l The difference between these two figures is subtle. In 
Fig. [25] the discriminant does not go below —2, and thus there are half as many gaps 
[96] . Direct observation of such a band structure would give unambiguous evidence of 
a PT-symmetric Hamiltonian [97] . Complex periodic potentials having more elaborate 
band structures have also been found [98] [991 HOOl HOlj - 

7.3. Quantum Brachistochrone 

We pointed out in Subsec. 16.41 that there is a similarity transformation that maps a 
non-Hermitian PT-symmetric Hamiltonian if to a Hermitian Hamiltonian h [see (11321) . 
The two Hamiltonians, H and h, have the same eigenvalues, but this does not mean 
they describe the same physics. To illustrate the differences between H and h, we 
show how to solve the quantum brachistochrone problem for PT-symmetric and for 
Hermitian quantum mechanics, and we show that the solution to this problem in these 
two formulations of quantum mechanics is not the same. 
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Figure 24. Discriminant A(£') plotted as a function of E for the real periodic potential 
V{x) — sin(a::). Although it cannot be seen in the figure, all local maxima lie above 
the line A — 2 and all local minima lie helow the line A = —2. The regions of 
energy E for which |A| < 2 correspond to allowed energy bands and the regions where 
|A| > 2 correspond to gaps. There are infinitely many gaps and these gaps become 
exponentially narrow as E increases. It is clear in this figure that the first maximum lies 
above 2. The first minimum occurs &t E = 2.313 8, where the discriminant has the value 
—2.003 878 7. The second maximum occurs aX E = 4.033 6, where the discriminant is 
2.000 007. Similar behavior is found for other potentials in the class V{x) = sin^^^^ {x). 



The fancy word brachistochrone means "shortest time." Thus, the quantum 
brachistochrone problem is defined as follows: Suppose we are given initial and 
final quantum states {ipi) and {ipp) in a Hilbert space. There exist infinitely many 
Hamiltonians H under which {ipj) evolves into \iPf) in some time t: 



The problem is to find the Hamiltonian H that minimizes the evolution time t subject 
to the constraint that uj, the difference between the largest and smallest eigenvalues of 
H, is held fixed. The shortest evolution time is denoted by r. 

In Hermitian quantum mechanics there is an unavoidable lower bound r on the 
time required to transform one state into another. Thus, the minimum time required 
to fiip unitarily a spin-up state into a spin-down state of an electron is an important 
physical quantity because it gives an upper bound on the speed of a quantum computer. 

In this paper we have shown that Hermitian quantum mechanics can be extended 
into the complex domain while retaining the reality of the energy eigenvalues, the 
unitarity of time evolution, and the probabilistic interpretation. It has recently been 
discovered that within this complex framework a spin-up state can be transformed 
arbitrarily quickly to a spin-down state by a simple non- Hermitian Hamiltonian [102] . 

Let us first show how to find the value of r for the case of a Hermitian Hamiltonian: 
This problem is easy because finding the optimal evolution time r requires only the 
solution to the simpler problem of finding the optimal evolution time for the 2x2 
matrix Hamiltonians acting on the two-dimensional subspace spanned by {ipi) and {ipp) 



iHt/h 



(147) 
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Figure 25. Discriminant A(iJ) plotted as a function of E for the complex periodic 
potential V{x) = «sin(x). For ii^ > 2 it is not possible to see any difference between 
this figure and Fig. [Ml However, although it cannot be seen, all local maxima in 
this figure lie above the line A = 2 and all local minima lie above the line A = — 2. 
This behavior is distinctly different from the behavior in Fig. [24l exhibited by the real 
periodic potential V{x) = sin a;. In stark contrast with real periodic potentials, for the 
potentials of the form V{x) — i sin^^^^(a;), while the maxima of the discriminant lie 
above +2, the minima of the discriminant lie above —2. Thus, for these potentials there 
are no antiperiodic wave functions. Lengthy and delicate numerical analysis verifies 
that for the potential V{x) — ism{x) the first three maxima of the discriminant A(£^) 
are located &t E ^ 3.966 428 4, E = 8.985 732 0, and E = 15.992 066 213 46. The 
value of the discriminant A(£;) at these energies is 2.000 007, 2.000 000 000 000 69, and 
2.000 000 000 000 000 000 04. The first two minima of the discriminant are located at 
E = 2.191 6 and 6.229 223 and at these energies 1\{E) has the values -1.995 338 6 
and —1.999 999 995 27. Similar behavior is found for the other complex periodic 
potentials in the class. 



|103] . To solve the Hermitian version of the two-dimensional quantum brachistochrone 
problem, we choose the basis so that the initial and final states are given by 

1^/) = J j and \^F)=(yl^. (148) 
where |ap + |6p = l. The most general 2x2 Hermitian Hamiltonian is 

H={ \, (149) 
\ re u I 

where the parameters r, s, u, and 6 are real. The eigenvalue constraint is 

uj^ = {s-uY + 4:r^. (150) 
The Hamiltonian H in f ll49p can be expressed in terms of the Pauli matrices 
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as H = |(s + u)l + icucr-n, where n = ^ (r cos ^, r sin ^, is a unit vector. The 
matrix identity 

exp(i(/) cr-n) = cos </) 1 + z sin </) cr-n (152) 
then simphfies the relation \iI)f) = e~^^'^^^\ipj) to 

/ ^os (i) - ^^ sin (i) \ 




^-i{s+u)t/{2h) 



(153) 



From the second component of (11531) we obtain |6| = ^ sin (|^), which gives the 



time required to transform the initial state: t = ^ arcsin (^). To optimize this relation 
over all r > 0, we note that fll50p implies that the maximum value of r is |a; and that 
at this maximum s = u. The optimal time is thus 

2h , , 

T = — arcsin . (154) 

u 

Note that if a = and 6 = 1 we have r = 27rh/uj for the smallest time required to 
transform (q) to the orthogonal state Cj*). 

For general a and b, at the optimal time r we have a = e^'^^I^^X — |6p and 
h = -ie^^/^lfele*^, which satisfies the condition |ap + |6p = 1 that the norm of the state 
does not change under unitary time evolution. The parameters s and 6 are determined 
by the phases of a and b. We set a = |a|e*'^''^^'^^ and b = |6|e*'^'"^(''^ and find that the 
optimal Hamiltonian is 

( ^arg(a) <£ j[arg{6)-arg(a)-f] \ 

2arcsm|b| 4 1 

c.; j[arg{6)-arg(a)-f] '^arg(a) I • V^'^^) 

4 2 arcsin |fe| / 

The overall phase of {ipp) is not physically relevant, so the quantity arg(a) may be chosen 
arbitrarily; we may thus assume that it is 0. We are free to choose arg(a) because there 
is no absolute energy in quantum mechanics; one can add a constant to the eigenvalues 
of the Hamiltonian without altering the physics. Equivalently, this means that the value 
of T cannot depend on the trace s + m of if. 

How do we interpret the result for r in (11541) ? While this equation resembles 
the time-energy uncertainty principle, it is really the much simpler statement that 
ratextime= distance. The constraint (11501) on if is a bound on the standard deviation 
AH of the Hamiltonian, where AH in a normalized state is given by [AH^ = 
{il)\H'^\ip) — {ip\H\ipy. The maximum value of AH is 0^/2. From the Anandan-Aharonov 
relation [104] . the speed of evolution of a quantum state is given by AH. The distance 
between the initial state \ipi) and the final state \ipp) is 6 = 2arccos(|('?/'i7'|'?/'7)|). Thus, 
the shortest time r for \ipj) to evolve to l^pp) = e~^^'^^^\il)i) is bounded below because the 
speed is bounded above while the distance is held fixed. The Hamiltonian H that realizes 
the shortest time evolution can be understood as follows: The standard deviation AH 
of the Hamiltonian in (I149P is r. Since AH is bounded by to maximize the speed 
of evolution (and minimize the time of evolution) we choose r = uj/2. 
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We now perform the same optimization for the complex 2x2 non-Hermitian VT 
Hamiltonian in fl93|) . Following the same procedure used for Hermitian Hamiltonians, 
we rewrite H in (1931) in the form H = (r cos^)l + ^u!(T-n, where n = ^(s, 0,2rsin6') is 
a unit vector and the squared difference between the energy eigenvalues [see (195!) is 

u'^ = is^ - ir^sin^ 9. (156) 

The condition of unbroken VT symmetry ensures the positivity of u;^. This equation 
emphasizes the key difference between Hermitian and PT-symmetric Hamiltonians: The 
corresponding equation (11501) for the Hermitian Hamiltonian has a sum of squares, while 
this equation for uj"^ has a difference of squares. Thus, Hermitian Hamiltonians exhibit 
elliptic behavior, which leads to a nonzero lower bound for the optimal time r. The 
hyperbolic nature of (I156P allows r to approach zero because, as we will see, the matrix 
elements of the PT-symmetric Hamiltonian can be made large without violating the 
energy constraint — E_ = oj. 

The PT-symmetric analog of the evolution equation (I153P is 




g— itr cosd/h 



COS a 



(157) 

V -^sin(g) ) 

We apply this result to the same pair of vectors examined in the Hermitian case: 
IV'/) = (q) and \'iI)f) = (i)- (Note that these vectors are not orthogonal with respect to 
the CVT inner product.) Equation (I157p shows that the evolution time to reach \iPf) 
from is t = (2a — 7r)h/uj. Optimizing this result over allowable values for a, we find 
that as a approaches ^vr the optimal time r tends to zero. 

Note that in the limit a — ^ ^vr we get cos a ^ 0. However, in terms of a, the energy 
constraint (11560 becomes tu^ = cos^ a. Since u is fixed, in order to have a approach 
^TT we must require that s ^ 1. It then follows from the relation sin a = (r/s) sin^ that 
|r| ~ \s\, so we must also require that r ^ 1. Evidently, in order to make r <^ 1, the 
matrix elements of the PT-symmetric Hamiltonian (1931) must be large. 

The result demonstrated here does not violate the uncertainty principle. Indeed, 
Hermitian and non-Hermitian PT-symmetric Hamiltonians share the properties that 
(i) the evolution time is given by 2'Kh/uj, and (ii) AH < uj/2. The key difference is that 
a pair of states such as (q) and are orthogonal in a Hermitian theory, but have 
separation 5 = vr — 2|a| in the PT-symmetric theory. This is because the Hilbert space 
metric of the PT-symmetric quantum theory depends on the Hamiltonian. Hence, it is 
possible to choose the parameter a to create a wormhole-like effect in the Hilbert space. 

A gedanken experiment to realize this effect in a laboratory might work as follows: A 
Stern-Gerlach filter creates a beam of spin-up electrons. The beam then passes through 
a 'black box' containing a device governed by a PT-symmetric Hamiltonian that flips 
the spins unitarily in a very short time. The outgoing beam then enters a second Stern- 
Gerlach device which verifies that the electrons are now in spin-down states. In effect, 
the black-box device is applying a magnetic field in the complex direction (s, 0, ir sin 6). 
If the field strength is sufficiently strong, then spins can be flipped unitarily in virtually 
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no time because the complex path joining these two states is arbitrary short without 
violating the energy constraint. The arbitrarily short alternative complex pathway 
from a spin- up state to a spin-down state, as illustrated by this thought experiment, is 
reminiscent of the short alternative distance between two widely separated space-time 
points as measured through a wormhole in general relativity. 

The results established here provide the possibility of performing experiments 
that distinguish between Hermitian and PT-symmetric Hamiltonians. If practical 
implementation of complex PT-symmetric Hamiltonians were feasible, then identifying 
the optimal unitary transformation would be particularly important in the design 
and implementation of fast quantum communication and computation algorithms. Of 
course, the wormhole-like effect we have discussed here can only be realized if it is 
possible to switch rapidly between Hermitian and PT-symmetric Hamiltonians by 
means of similarity transformations. It is conceivable that so much quantum noise 
would be generated that there is a sort of quantum protection mechanism that places 
a lower bound on the time required to switch Hilbert spaces. If so, this would limit the 
applicability of a Hilbert-space wormhole to improve quantum algorithms. 

7.4- Supersymmetric VT -Symmetric Hamiltonians 

After the discovery of PT-symmetric Hamiltonians in quantum mechanics, it was 
proposed that VT symmetry might be combined with supersymmetry [105j in the 
context of quantum field theory. In Ref. [105] it was shown that one can easily construct 
two-dimensional supersjTiimetric quantum field theories by introducing a PT-symmetric 
superpotential of the form 5(0) = —ig{i(f)Y^'^ . The resulting quantum field theories 
exhibit a broken parity symmetry for all 5 > 0. However, supersymmetry remains 
unbroken, which is verified by showing that the ground-state energy density vanishes 
and that the fermion-boson mass ratio is unity. Many papers have subsequently worked 
on the subject of "PT-symmetric supersymmetric quantum mechanics. (See Ref. |106j .) 
A particularly interesting paper by Dorey et al. examines the connection between 
supersymmetry and broken VT symmetry in quantum mechanics jl07j . 

7. 5. Other Quantum- Mechanical Applications 

There are many additional quantum-mechanical applications of non-Hermitian VT- 
invariant Hamiltonians. In condensed matter physics Hamiltonians rendered non- 
Hermitian by an imaginary external field have been introduced to study delocalization 
transitions in condensed matter systems such as vortex flux-line depinning in type- 
II superconductors [108] . In this Hatano- Nelson model there is a critical value of the 
anisotropy (non-Hermiticity) parameter below which all eigenvalues are real [109j . In the 
theory of Reaction-Diffusion systems, many models have been constructed for systems 
described by matrices that can be non-Hermitian |110l 1111] and with the appropriate 
definition of the V and T operators, these systems can be shown to be VT symmetric. 
Finally, we mention a recent paper by Hibberd et al. who found a transformation that 
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maps a Hamiltonian describing coherent superpositions of Cooper pairs and condensed 
molecular bosons to one that is PT-symmetric |112j . 

8. PT-Symmetric Quantum Field Theory 

Quantum-mechanical theories have only a finite number of degrees of freedom. Most 
of the PT-symmetric quantum-mechanical models discussed so far in this paper have 
just one degree of freedom; that is, the Hamiltonians for these theories are constructed 
from just one pair of dynamical variables, x and p. In a quantum field theory the 
operators x{t) and p{t) are replaced by the quantum fields !f{'x,t) and 7r(x, t), which 
represent a continuously infinite number of degrees of freedom, one for each value of the 
spatial variable x. Such theories are vastly more complicated than quantum-mechanical 
theories, but constructing quantum field theories that are non-Hermitian and VT 
symmetric is straightforward. For example, the quantum-field-theoretic Hamiltonians 
that are analogous to the quantum-mechanical Hamiltonians in ( ITOj) and (fTTj) are 

n = i7r2(x, t) + i[Vx¥^(x, t)]^ + |/iV'(x, t) + ig^\^, t) (158) 

and 

n = |7r2(x,t) + \[VM^,t)f + |/"V'(x,t) - \gv'M. (159) 

As in quantum mechanics, where the operators x and p change sign under parity 
refiection P, we assume that the fields in these Hamiltonians are pseudoscalars so that 
they also change sign: 

P¥?(x, t)V = -y?(-x, t), P7r(x, t)V = -7r(-x, t). (160) 

Quantum field theories like these that possess VT symmetry exhibit a rich variety of 
behaviors. Cubic field-theory models like that in (1158P are of interest because they arise 
in the study of the Lee- Yang edge singularity [HI [821 [83] and in Reggeon field theory 
[85| [86] . In these papers it was asserted that an i(p^ field theory is nonunitary. However, 
by constructing the C operator, we argue in Subsec. 18. II that this quantum field theory 
is, in fact, unitary. We show in Subsec. 18.31 how VT symmetry eliminates the ghosts 
in the Lee Model, another cubic quantum field theory. The field theory described by 
(11591) is striking because it is asymptotically free, as explained in Subsec. 18. 4[ We also 
examine PT-symmetric quantum electrodynamics in Subsec. 18.51 the "PT-symmetric 
Thirring and Sine-Gordon models in Subsec. 18.61 and gravitational and cosmological 
theories in Subsec. 18. 7[ Last, we look briefiy at PT-symmetric classical field theories in 
Subsec. EB 

8.1. iif^ Quantum Field Theory 

In courses on quantum field theory, a scalar gifr" theory is used as a pedagogical example 
of perturbative renormalization even though this model is not physically realistic because 
the energy is not bounded below. However, by calculating C perturbatively for the case 
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when the couphng constant g = it is pure imaginary, one obtains a fully acceptable 
Lorentz invariant quantum field theory. This calculation shows that it is possible to 
construct perturbatively the Hilbert space in which the Hamiltonian for this cubic scalar 
field theory in [D + l)-dimensional Minkowski space-time is self-adjoint. Consequently, 
such theories have positive spectra and exhibit unitary time evolution. 

In this subsection we explain how to calculate perturbatively the C operator 
for the quantum- field-theoretic Hamiltonian in (11581) [71] . We apply the powerful 
algebraic techniques explained in Subsecs. 16.11 and 16.21 for the calculation of the C 
operator in quantum mechanics. As in quantum mechanics, we express C in the form 
C = exp {eQi + e^Qs + . . .)V, where now Q2n+i = 0, 1, 2, . . .) are real functionals 
of the fields (^x and tt^. To find Q„ for Ti in (11581) we must solve a system of operator 
equations. 

We begin by making an ansatz for Qi analogous to the ansatz used in (11201) : 

Qi = j j j dyidydz [Mi^^y^-^n^-Ky-K^ + iVx(yz)V^y7rx</2z) • (161) 

In quantum mechanics M and are constants, but in field theory they are functions. 
The notation M(xyz) indicates that this function is totally symmetric in its three 
arguments, and the notation Ax(yz) indicates that this function is symmetric under 
the interchange of the second and third arguments. The unknown functions M and 
are form factors; they describe the spatial distribution of three-point interactions of 
the field variables in Qi. The nonlocal spatial interaction of the fields is an intrinsic 
property of C. (Note that we have suppressed the time variable t in the fields and that 
we use subscripts to indicate the spatial dependence.) 

To determine M and A^ we substitute Qi into the first equation in (11191) . namely 
[Hq, Qi] = —2Hi, which now takes the form 

j d^7il + j I rfxrfy^xG'x>y,Qi = -42 j d^^l, (162) 

where the inverse Green's function is given by G^^ = (/i^ - V^)5(x - y). We obtain the 
following coupled system of partial differential equations: 

(/i^ - V^)Arx(yz) + (/i^ - Vj)Ary(xz) + (/i^ - V^)Ar,(xy) 
= -65(x-y)5(x-z), 

A^x(yz) + Ary(xz) = 3(/i2 - V^)M(wxy). (163) 

We solve these equations by Fourier transforming to momentum space and get 

4 f f g«(x-y)-p+«(x-z)-q 
= / / "Pi'i ^(p_q) ■ (164) 

where -D(p, q) = 4[p^q^ — (p ■ q)^] + 4/i^(p^ + P ■ Q + q^) + is positive, and 

A^x(yz) = - 3 (Vy ■ V. + i/i') M(xyz). (165) 

For the special case of a (1 + l)-dimensional quantum field theory the integral in (11651) 
evaluates to M(xyz) = — Ko(/i-R)/ (a/Stt/x^), where Kq is the associated Bessel function 
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and = |[(x — y)^ + (y — z)^ + (z — x)^]. This completes the calculation of C to first 
order in perturbation theory. 

We mention finally that the C operator for this cubic quantum field theory 
transforms as a scalar under the action of the homogeneous Lorentz group |113j . In 
Ref. |113j it was argued that because the Hamiltonian Hq for the unperturbed theory 
{g = 0) commutes with the parity operator V, the intrinsic parity operator Vi in the 
noninteracting theory transforms as a Lorentz scalar. (The intrinsic parity operator Vi 
and the parity operator V have the same effect on the fields, except that Vi does not 
reverse the sign of the spatial argument of the field. In quantum mechanics V and Vi 
are indistinguishable.) When the coupling constant g is nonzero, the parity symmetry 
of H is broken and Vi is no longer a scalar. However, C is a. scalar. Since lim^^^o C = "Pi, 
one can interpret the C operator in quantum field theory as the complex extension of the 
intrinsic parity operator when the imaginary coupling constant is turned on. This means 
that C is frame-invariant and it shows that the C operator plays a truly fundamental 
role in non-Hermitian quantum field theory. 



8.2. Other Quantum Field Theories Having Cubic Interactions 

We can repeat the calculations done in Subsec. 18.11 for cubic quantum field theories 
having several interacting scalar fields [71], [72] . For example, consider the case of two 
scalar fields if^^ and (p^^ whose interaction is governed by 

H = i/i^) + //f +^e|d^ [^^'f , (166) 

which is the analog of the quantum- mechanical theory described by H in fll23p . Here, 

4'^ = i 1 rfx [n^J^f + ljjd^dy [Gg]-V^^V?'^- (167) 

To determine C to order e we introduce the ansatz 

r(l) ^^(l),^(l) _^ ,^{1)^(1)^1 ,^(2) 

(168) 



(2) ,^{2),J1),J1) , M , .^(2)^(1)^(1) 



y ^z 



where Mx(yz), -^x{yz)' ^'^'^ x(yz) ^'^^ unknown functions and the parentheses indicate 
symmetrization. We get 

^x(yz) = — GmiRl) G,_i2{R2), 

<yz) = - 5(2x-y-z)G™(i?i), 

^Syz) = i^(2x - y - z)GmiRi) - 5(y - z)G^,(i?2), (169) 

where G^{r) = j^{^)^^'^K_ij^£)/2ifJ'r) with r = |r| is the Green's function in D- 
dimensional space, m? = — jul, Rf = {y — z)^, and = j{2x — y — z)^. 
For three interacting scalar fields whose dynamics is described by 
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which is the analog of H in (11251) . we make the ansatz 



(171) 



The solutions for the unknown functions are as follows: Mxyz is given by the integral 
(11641) with the more general formula D{p, q) = 4[p^q^ — (p ■ q)^] + 4[yU^(q^ + P " q) + 
/i2(p2 + p . q) _ . qj ^ ^4 ^-^^Yi /i^ = 2/i2/i2 + 2/^2^2 ^ 2ij,lnl - nf- The 

coefficients are expressed as derivatives acting on M: 



Once again, the calculation of C shows that these cubic field theories are fully consistent 
quantum theories lffl 

8.3. The Lee model 

In this subsection we will show how to use the tools that we have developed to study non- 
Hermitian "PT-symmetric quantum theories to examine a famous model quantum field 
theory known as the Lee model. In 1954 the Lee model was proposed as a quantum field 
theory in which mass, wave function, and charge renormalization could be performed 
exactly and in closed form [H I114[ I115i 1116] . We discuss the Lee model here because 
when the renormalized coupling constant is taken to be larger than a critical value, 
the Hamiltonian becomes non-Hermitian and a (negative- norm) ghost state appears. 
The appearance of the ghost state was assumed to be a fundamental defect of the Lee 
model. However, we show that the non-Hermitian Lee-model Hamiltonian is actually 
VT symmetric. When the states of this model are examined using the C operator, the 
ghost state is found to be an ordinary physical state having positive norm |117] . 

The idea for studying the Lee model as a non-Hermitian Hamiltonian is due to 
Kleefeld, who was the first to point out this transition to VT symmetry |118] . His work 
gives a comprehensive history of non-Hermitian Hamiltonians. 

The Lee model has a cubic interaction term that describes the interaction of three 
spinless particles called V, N, and 6. The V and particles are fermions and behave 
roughly like nucleons, and the 6 particle is a boson and behaves roughly like a pion. In 
the model a V may emit a 6, but when it does so it becomes an A^: V ^ N + 6. Also, 
an A^ may absorb a 6, but when it does so it becomes aV: N + 9 ^ V . 

The solvability of the Lee model is based on the fact that there is no crossing 
symmetry. That is, the A^ is forbidden to emit an anti-^ and become a V. Eliminating 

tfln Refs. |113l 1117] it is shown that the C operator in quantum field theory has the form C = e'^Vi, 
where Vi is the intrinsic parity reflection operator. The difference between V and Vi is that Vi does 
not reflect the spatial arguments of the flelds. For a cubic interaction Hamiltonian this distinction is 
technical. It does not affect the final result for the Q operator in (|16ip . (|168p . and (|17ip . 




(172) 
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crossing symmetry makes the Lee model solvable because it introduces two conservation 
laws. First, the number of N quanta plus the number of V quanta is fixed. Second, 
the number of N quanta minus the number of 9 quanta is fixed. These two highly 
constraining conservation laws decompose the Hilbert space of states into an infinite 
number of noninteracting sectors. The simplest sector is the vacuum sector. Because 
of the conservation laws, there are no vacuum graphs and the bare vacuum is the 
physical vacuum. The next two sectors are the one-6'-particle and the one-A^-particle 
sector. These two sectors are also trivial because the two conservation laws prevent any 
dynamical processes from occurring there. As a result, the masses of the N particle and 
of the 9 particle are not renormalized; that is, the physical masses of these particles are 
the same as their bare masses. 

The lowest nontrivial sector is the V/N9 sector. The physical states in this sector 
of the Hilbert space are linear combinations of the bare V and the bare N9 states, 
and these states consist of the one-physical- V-particle state and the physical A^-^- 
scattering states. To find these states one can look for the poles and cuts of the 
Green's functions. The renormalization in this sector is easy to perform. Following 
the conventional renormalization procedure, one finds that the mass of the V particle is 
renormalized; that is, the mass of the physical V particle is different from its bare mass. 
In the Lee model one calculates the unrenormalized coupling constant as a function of 
the renormalized coupling constant in closed form. There are many ways to define the 
renormalized coupling constant. For example, in an actual scattering experiment one 
could define the square of the renormalized coupling constant g"^ as the value of the A^^ 
scattering amplitude at threshold. 

The intriguing aspect of the Lee model is the appearance of a ghost state in the 
V/N9 sector. This state appears when one performs coupling-constant renormalization. 
Expressing gfg, the square of the unrenormalized coupling constant, in terms of g"^, 
the square of the renormalized coupling constant, one obtains the graph in Fig. [H 
In principle, the is a physical parameter whose numerical value is determined by a 
laboratory experiment. If g"^ is measured to be near 0, then from Fig. [26] we see that gQ 
is also small. However, if the experimental value of g'^ is larger than this critical value, 
then the square of the unrenormalized coupling constant is negative. In this regime go 
is imaginary and the Hamiltonian is non-Hermitian. Moreover, a new state appears in 
the V/N9 sector, and because its norm is negative, the state is called a ghost. Ghost 
states are unacceptable in quantum theory because their presence signals a violation of 
unitarity and makes a probabilistic interpretation impossible. 

There have been many unsuccessful attempts to make sense of the Lee model as 
a physical quantum theory in the ghost regime [H I115[ 1116] . However, the methods of 
PT-symmetric quantum theory enable us to give a physical interpretation for the V/N9 
sector of the Lee model when go becomes imaginary and H becomes non-Hermitian. The 
Lee model is a cubic interaction and we have already shown in Subsecs. 18. II and ] 8.2l how 
to make sense of a Hamiltonian in which there is a cubic interaction multiplied by an 
imaginary coupling constant. The procedure is to calculate the C operator and to use 
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it to define a new inner product when the Hamiltonian is non-Hermitian. 

For simphcity, we focus here on the quantum-mechanical Lee model; the results for 
the field-theoretic Lee model in Ref. [11 7j are qualitatively identical. The Hamiltonian 
for the quantum-mechanical Lee model is 

+ mea^a+ (V'^Na + o^N'^V^ -(173) 

The bare states are the eigenstates of Hq and the physical states are the eigenstates of 
the full Hamiltonian H. The mass parameters mj^ and uiq represent the physical masses 
of the one-iV-particle and one-^-particle states because these states do not undergo mass 
renormalization. However, my^ is the bare mass of the V particle. 

We treat the V, N, and 9 particles as pseudoscalars. To understand why, recall 
that in quantum mechanics the position operator x = {a + a'^) / \/2 and the momentum 
operator p = i{a^ — a)/^/2 both change sign under parity reflection V: 

VxV = -X, VpV = -p. (174) 

Thus, VVV = -V, VNV = -N, VaV = -a, VV^V = -V^, VN^V = 
Va'^V = —a^. Under time reversal T, p and i change sign but x does not: 

TpT = —p, TiT = —i, TxT = x. (175) 

Thus, rVT = V, TNT = N, TaT = a, TV^T = V\ TN^T = N\ Ta^T = al 

When the bare coupling constant go is real, H in ( 1173^ is Hermitian: = H. 

When go is imaginary, go = iXo (Aq real), H is not Hermitian, but by virtue of the 

above transformation properties, H is PT-symmetric: H^^ = H. 

We assume first that go is real so that H is Hermitian and we examine the simplest 

nontrivial sector of the quantum- mechanical Lee model; namely, the V/N6 sector. 




Figure 26. Square of the unrenormalized coupling constant, g^, plotted as a function 
of the square of the rcnormalized coupling constant g^. Note that g^ — when = 0, 
and as increases from so does 5q. However, as increases past a critical value, 
gQ abruptly becomes negative. In this regime go is imaginary and the Hamiltonian is 
non-Hermitian. 
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We look for the eigenstates of the Hamihonian H in the form of hnear combinations 
of the bare one- ^-particle and the bare one-A^-one-^-particle states. There are two 
eigenf unctions. We interpret the eigenfunction corresponding to the lower-energy 
eigenvalue as the physical one-\^-particle state, and we interpret the eigenfunction 
corresponding with the higher-energy eigenvalue as the physical one-iV-one-^-particle 
state. (In the field-theoretic Lee model this higher-energy state corresponds to the 
continuum of physical N-9 scattering states.) Thus, we make the ansatz 

|y) = cn|l,0,0) + ci2|0,l,l), |iV^^) =C2i|l,0,0) + C22|0,l,l), (176) 

and demand that these states be eigenstates of H with eigenvalues my (the renormalized 
l^-particle mass) and E^e- The eigenvalue problem reduces to a pair of elementary 
algebraic equations: 

ciimvo + 0125^0 = ciimy, CaiS^o + C22 ("^tv + me) = C22^7ve- (177) 
The solutions to (11771) are 



my =^ [rriN + mg + mvo - \/ A^o + 



En9 = ^[mN + mg + mvo + \/ fJ^l + %o ) > (178) 

where /io = rriN + mo — my^- Notice that my, the mass of the physical V particle, is 
different from myp, the mass of the bare V particle, because the V particle undergoes 
mass renormalization. 

Next, we perform wave-function renormalization. Following Barton we define 
the wave-function renormalization constant Zy by ^/Zy = {0\V\V). This gives 



Zy' = \go^\['^^+H - /"o j . (179) 

Finally, we perform coupling-constant renormalization. Again, following Barton we 
note that ^/Zy is the ratio between the renormalized coupling constant g and the bare 
coupling constant g^ [Ij. Thus, g'^ I g^ = Zy. Elementary algebra gives the bare coupling 
constant in terms of the renormalized mass and coupling constant: 

9l = 9V{l-9Vf^'), (180) 

where /i is defined as = mN + mg — my. We cannot freely choose g because the 
value of g is, in principle, taken from experimental data. Once g has been determined 
experimentally, we can use ( l2Ti) to determine go. The relation in ( l2Ti) is plotted in 
Fig. [26l This figure reveals a surprising property of the Lee model: If g is larger than 
the critical value fi, then the square of 5^0 is negative and 5^0 is imaginary. 

As g approaches its critical value from below, the two energy eigenvalues in (11781) 
vary accordingly. The energy eigenvalues are the two zeros of the secular determinant 
f{E) obtained from applying Cramer's rule to (11771) . As g (and ^^o) increase, the energy 
of the physical N6 state increases. The energy of the A^^ state becomes infinite as 
g reaches its critical value. As g increases past its critical value, the upper energy 
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eigenvalue goes around the bend; it abruptly jumps from being large and positive to 
being large and negative. Then, as g continues to increase, this energy eigenvalue 
approaches the energy of the physical V particle from below. 

When g increases past its critical value, the Hamiltonian H in (11731) becomes non- 
Hermitian, but its eigenvalues in the V/N9 sector remain real. (The eigenvalues remain 
real because H becomes VT symmetric. All cubic PT-symmetric Hamiltonians that 
we have studied have been shown to have real spectra.) However, in the "PT-symmetric 
regime it is no longer appropriate to interpret the lower eigenvalue as the energy of the 
physical NO state. Rather, it is the energy of a new kind of state \G) called a ghost. As 
is shown in Refs. |115l 11161 H] , the Hermitian norm of this state is negative. 

A physical interpretation of the ghost state emerges easily when we use the 
procedure developed in Ref. [117] . We begin by verifying that in the PT-symmetric 
regime, where go is imaginary, the states of the Hamiltonian are eigenstates of the VT 
operator, and we then choose the multiplicative phases of these states so that their VT 
eigenvalues are unity: 

VT\G) = \G), VT\V) = \V). (181) 

It is then straightforward to verify that the VT norm of the V state is positive, while 
the VT norm of the ghost state is negative. 

We now follow the procedures described in Sec. [6] to calculate C. We express 
the C operator as an exponential of a function Q multiplying the parity operator: 
C = exp [QiV'^ ,V; N'^ , N; a\ a)~\ V. We then impose the operator equations = 1, 
[C, VT] = 0, and [C, H] = 0. The condition = 1 gives 

Q{V\ V] N\ N; a\ a) = -Q{-V\ -V; -N\ -N; -a\ -a). (182) 

Thus, QiV^ , V; N'^ , N; , a) is an odd function in total powers of , V, N'^ , N, , and 
a. Next, we impose the condition [C,VT] = and obtain 

Q{V\ V; N\ N; a^ a) = Q*i-V\ -V; -N\ -N; -a^ -a), (183) 

where * denotes complex conjugation. 

Last, we impose the condition that C commutes with H, which requires that 

[e«,/7o] =go[e'^,H,]^. (184) 

Although in Subsecs. 18.11 and 18.21 we were only able to find the C operator to leading 
order in perturbation theory, for the Lee model one can calculate C exactly and in closed 
form. To do so, we seek a solution for Q as a formal Taylor series in powers of go: 

oo 

Q = J29o'"''Q2n+i. (185) 

n=0 

Only odd powers of go appear in this series, and Q2n+i are all ant i- Hermitian: Qln+i = 
-Q2n+i- From (11851) we get 

22n+l 

g^n+i = (-l)" (2n+l)/xg"+^ {V^Nan^ - n^a^N^V) , (186) 
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where hq = a^a is the number operator for ^-particle quanta. 

We then sum over all Q2n+i and obtain the exact result that 

1 f2gQJn^\ 1 {2goJn^ 



Q = V^Nc 



larctan 



: arctan 



i^N^V. (187) 



'ne V /io y ^/ne \ jiQ 

We exponentiate this result to obtain the C operator. The exponential of Q simplifies 
considerably because we are treating the V and N particles as fermions and therefore 
we can use the identity nyj^ = nv,N- The exact result for e'^ is 



1 



/ionAr(l-ny) i2onv{l-nN) 



+ 



+ 4c/o^e ^/J4T^gf(j^^TT) 



ny — + nyUN 



We can also express the parity operator V in terms of number operators: 



p = ^iAnv+n^+ne) = {I _ 2ny) (1 - 2nN} 6 

Combining e'^ and P, we obtain the exact expression for C: 



189) 



C 



1 



2go 



2^70 



^tjyt^ _^ /iprajv (1 - ny) _^ /ipny (1 - rijv) 



— ny — + nyn^ 



;i - 2ny)(l - 2njv) e*^'^^ (190) 



Using this C operator to calculate the CVT norm of the V state and of the 
ghost state, we find that these norms are both positive. Furthermore, the the time 
evolution is unitary. This establishes that with the proper definition of the inner 
product the quantum-mechanical Lee model is a physically acceptable and consistent 
quantum theory, even in the ghost regime where the unrenormalized coupling constant 
is imaginary. The procedure of redefining the inner product to show that the ghost state 
is a physical state is a powerful technique that has been used by Curtright et al. and 
by Ivanov and Smilga for more advanced problems |119[ 1120] . 



8.4- The Higgs Sector of the Standard Model of Particle Physics 

The distinguishing features of the —g^'^ quantum field theory in fll59p are that its 
spectrum is real and bounded below, it is perturbatively renormalizable, it has a 
dimensionless coupling constant in four-dimensional space-time, and it is asymptotically 
free [121J. The property of asymptotic freedom was established many years ago by 
Symanzik [122J, as has been emphasized in a recent paper by Kleefeld |123] . As 
Kleefeld explains, a +g(p^ theory in four- dimensional space-time is trivial because it 
is not asymptotically free. However, Symanzik proposed a "precarious" theory with a 
negative quartic coupling constant as a candidate for an asymptotically free theory of 
strong interactions. Symanzik used the term "precarious" because the negative sign of 
the coupling constant suggests that this theory is energetically unstable. However, as we 
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have argued in this paper, imposing PT-symmetric boundary conditions (in this case 
on the functional-integral representation of the quantum field theory) gives a spectrum 
that is bounded below. Thus, Symansik's proposal of a nontrivial theory is resurrected. 

The —g^p'^ quantum field theory has another remarkable property. Although the 
theory seems to be parity invariant, the PT-symmetric boundary conditions violate 
parity invariance, as explained in Subsec. 12. 7[ Hence, the one-point Green's function 
(the expectation value of the field ip) does not vanish. (Techniques for calculating this 
expectation value are explained in [124] .) Thus, a nonzero vacuum expectation value 
can be achieved without having to have spontaneous symmetry breaking because parity 
symmetry is permanently broken. These properties suggest that a —g^p'^ quantum field 
theory might be useful in describing the Higgs sector of the standard model. 

Perhaps the Higgs particle state is a consequence of the field-theoretic parity 
anomaly in the same way that the quantum-mechanical parity anomaly described in 
Subsecs. 12.71 and 12.81 gives rise to bound states. Recent research in this area has been 
done by Jones et ai, who studied transformations of functional integrals |125] . and by 
Meisinger and Ogilvie, who worked on large- iV approximations and matrix models [126j . 

8.5. VT -Symmetric Quantum Electrodynamics 

If the unrenormalized electric charge e in the Hamiltonian for quantum electrodynamics 
were imaginary, then the Hamiltonian would be non-Hermitian. However, if one specifies 
that the potential in this theory transforms as an axial vector instead of a vector, 
then the Hamiltonian becomes VT symmetric [127] . Specifically, we assume that the 
four-vector potential and the electromagnetic fields transform under V like 

V : E ^ E, B ^ -B, A ^ A, A" (191) 

Under time reversal, the transformations are assumed to be conventional: 

r : E ^ E, B ^ -B, A ^ -A, A° (192) 

The Lagrangian of the theory possesses an imaginary coupling constant in order that it 
be invariant under the product of these two symmetries: 

C = -iFf'^F^^ + ^tph^Y^d^^ + ^mtp^j^ip + ie^^j^Y^jAf,. (193) 

The corresponding Hamiltonian density is then 

n = \{E^ + B^) + [7%'= {-iVk + leAk) + m7°] ^. (194) 

The Lorentz transformation properties of the fermions are unchanged from the usual 
ones. Thus, the electric current appearing in the Lagrangian and Hamiltonian densities, 
j/* = ip^'j^^'yt^ilj^ transforms conventionally under both V and T: 

Vf{^,t)V= j (-x,t), Tf{^,t)T= j (x,-t). (195) 

Because its interaction is cubic, this non-Hermitian theory of "electrodynamics" 
is the analog of the spinless iip^ quantum field theory discussed in Subsec. 18. 1[ VT- 
symmetric electrodynamics is especially interesting because it is an asymptotically free 
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theory (unlike ordinary electrodynamics) and because the sign of the Casimir force is 
the opposite of that in ordinary electrodynamics |127[ 1128] . This theory is remarkable 
because finiteness conditions enable it to determine its own coupling constant [128] . 

The C operator for PT-symmetric quantum electrodynamics has been constructed 
perturbatively to first order in e [129] . This construction is too technical to describe 
here, but it demonstrates that non-Hermitian quantum electrodynamics is a viable 
and consistent unitary quantum field theory. PT-symmetric quantum electrodynamics 
is more interesting than an quantum field theory because it possesses many 
of the features of conventional quantum electrodynamics, including Abelian gauge 
invariance. The only asymptotically free quantum field theories described by Hermitian 
Hamiltonians are those that possess a non- Abelian gauge invariance; VT symmetry 
allows for new kinds of asymptotically free theories, such as the —ip'^ theory discussed 
in Subsec. 18.41 that do not possess a non- Abelian gauge invariance. 

8.6. Dual VT- Symmetric Quantum Field Theories 

Until now we have focused on bosonic "PT-symmetric Hamiltonians, but it is just as 
easy to construct fermionic PT-symmetric Hamiltonians. We look first at free theories. 
The Lagrangian density for a conventional Hermitian free fermion field theory is 

C{x,t) =^{x,t){i0 -m)tp{x,t) (196) 

and the corresponding Hamiltonian density is 

7^(x, t) = ^(x, t){-if+ m)V'(x, t), (197) 

where %l){ii.,t) = ip^{-iL,t)'jo. 

In (1 + l)-dimensional space-time we adopt the following conventions: 70 = o"i 
and 7i = ia2, where the Pauli a matrices are given in f llSip . With these definitions, 
we have 7q = 1 and •yf = —1. We also define 75 = 7071 = 0-3, so that 7I = 1. The 
parity-reflection operator V has the effect 

V^ix, t)V = 7o^(-x, t), Vi){x, t)V = ^(-x, t)7o. (198) 

The effect of the time-reversal operator T, 

rv'(x,t)r = 7o^(x,-t), r^(x,t)r = V'(x,-t)7o, (i99) 

is similar to that of P, except that T is antilinear and takes the complex conjugate of 
complex numbers. From these definitions the Hamiltonian H = J dxT-C{x,t), where Ti 
is given in (11971) . is Hermitian: H = W . Also, H is separately invariant under parity 
reflection and under time reversal: VHV = H and THT = H. 

We can construct a non-Hermitian fermionic Hamiltonian by adding a 75-dependent 
mass term to the Hamiltonian density in (I197P : 

n{x, t) = tpix, t){-i'f+ mi + m275)^/'(x, t) (ma real). (200) 

The Hamiltonian H = J dxT-C{x,t) associated with this Hamiltonian density is not 
Hermitian because the m2 term changes sign under Hermitian conjugation. This sign 
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change occurs because 70 and 75 anticommute. Also, H is not invariant under V or 
under T separately because the m2 term changes sign under each of these reflections. 
However, H is invariant under combined V and T reflection. Thus, H is PT-symmetric. 

To see whether the VT symmetry of H is broken or unbroken, we must check to 
see whether the spectrum of H is real. We do so by solving the field equations. The 
field equation associated with Ti. in (12001) is 

(i^-mi -m275) V^(x,t) = 0. (201) 

If we iterate this equation and use 0^ = 9"^, we obtain 

{d^ + fi^)^{x,t) = 0, (202) 

which is the two-dimensional Klein-Gordon equation with /i^ = ml — m\. The physical 
mass that propagates under this equation is real when the inequality 

m\ > ml (203) 

is satisfied. This condition defines the two-dimensional parametric region of unbroken 
VT symmetry. When (12031) is not satisfied, VT symmetry is broken. At the boundary 
between the regions of broken and unbroken VT symmetry (the line m2 = 0), the 
Hamiltonian is Hermitian. [The same is true in quantum mechanics. Recall that for the 
Hamiltonian in (fT2l) the region of broken (unbroken) VT symmetry is e < (e > 0). At 
the boundary e = of these two regions the Hamiltonian is Hermitian.] 

The C operator associated with the PT-symmetric Hamiltonian density Ti. in (12001) 
IS given by C = e^V, where pH] 



Q= —tanh J dxilj{x,t)'yiip{x,t) 

= -tanh-^^y dxip^{x,t)-f5ip{x,t). (204) 

The inverse hyperbolic tangent in this equation requires that \e\ < 1, or equivalently 
that mf > m\, which corresponds to the unbroken region of VT symmetry. We use 
(12041) to construct the equivalent Hermitian Hamiltonian h as in (I132p : 



h = exp |tanh J dxijj'^{x,t)'j5ip{x,t) 



H 



X exp 



ifanh ^£ I dxip'^ {x,t)'^^il){x,t) 



(205) 



The commutation relations [75,70] = —271 and [75,71] = —270 simplify h in (I205p : 

h= dxij{x,t){-tf+fi)'^{x,t), (206) 



where /x^ = m^(l — e"^) = m\ — m^, in agreement with (12021) . Replacing H hj h changes 
the 75-dependent mass term mtp{l + £75)^/' fo a conventional fermion mass term fiipip. 
Thus, the non-Hermitian PT-symmetric Hamiltonian density in ( 1200p is equivalent to 
the Hermitian Hamiltonian density in (I197p with m replaced by fi. 
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If we introduce a four-point fermion interaction term in (11961] . we obtain the 
Lagrangian density for the massive Thirring model in (1 + 1) dimensions: 

C = ^{i$-m)ij + \g{ijYi^){i^li.'^). (207) 
whose corresponding Hamiltonian density is 

n = m)tp - (208) 

We can then construct the PT-symmetric Thirring model 

n = ij{-if+ m + em-f5)^ - |^(V'7^V)(V'7mV') (209) 

by introducing a 75-dependent mass term. The additional term is non-Hermitian 
but PT-symmetric because it is odd under both parity reflection and time reversal. 
Remarkably, the Q operator for the interacting case g ^ is identical to the Q 
operator for the case g = because in (1 + l)-dimensional space the interaction term 
[ijj'y'^ijj){ilj'j^ilj) commutes with the Q in (12041) [130J. Thus, the non-Hermitian VT- 
symmetric Hamiltonian density in (12091) is equivalent to the Hermitian Hamiltonian 
density in (I208p with the mass m replaced by /i, where //^ = m^(l — e'^) = ml — m\. 

The same holds true for the (3 + l)-dimensional interacting Thirring model by 
virtue of the commutation relation [75,707/^] = 0, but because this higher-dimensional 
field theory is nonrenormalizable, the Q operator may only have a formal significance. 

In (1 + 1) dimensions the massive Thirring Model ( ]207p is dual to the (1 + 1)- 
dimensional Sine-Gordon model |131] . whose Lagrangian density is 

£= i(p^y,)2 + ^2^-2(cosA(/?-l), (210) 

whose corresponding Hamiltonian density is 

H = ivr^ + \{V^f + m^\~'^{l - cos \^), (211) 

where 7r(x,t) = dQip{x,t) and where in (1 + l)-dimensional space V'^i^x.t) is just 
Piip{x,t). The duality between the Thirring model and the Sine-Gordon model is 
expressed as an algebraic relationship between the coupling constants g and A: 

- = (212) 

Note that the free fermion theory [g = 0) is equivalent to the Sine-Gordon model with 
the special value for the coupling constant = Att. 

The PT-symmetric extension (12091) of the modified Thirring model is, by the same 
analysis, dual to a modified Sine-Gordon model with the Hamiltonian density 

H = ivr^ + i(V(/?)2 + m^X-\l - cos Xif - ie sin Xif), (213) 

which is PT-symmetric and not Hermitian. The Q operator for this Hamiltonian is 

Q = ^tsinh~^e [dxiT{x,t). (214) 



(215) 



A 

Thus, the equivalent Hermitian Hamiltonian h is given by 



h = exp 



— tanh ^ejdxn{x,t) H exp — tanh ^ejdx7c{x,t) 
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Note that the operation that transforms H to h has the same effect as shifting the boson 
field by an imaginary constant: 

(p ^ ip + ^t&nh-^e. (216) 
A 

Under this transformation the interaction term m^A~^(l — cosXip — zesinAy?) in (12131) 
becomes — m^A~^(l — e^) cos Ay?, apart from an additive constant. Hence, h is the 
Hamiltonian for the Hermitian Sine-Gordon model, but with mass /i given by /i^ = 
m^(l — e^) = ml — 1712- This change in the mass is the same as for the Thirring model. 
Being Hermitian, h is even in the parameter e that breaks the Hermiticity of H. 

The idea of generating a non-Hermitian but PT-symmetric Hamiltonian from a 
Hermitian Hamiltonian by shifting the field operator as in (12161) . first introduced in 
the context of quantum mechanics in Ref. |132j . suggests a way to construct solvable 
fermionic PT-invariant models whenever there is a boson-fermion duality. 

8.7. VT -Symmetric Gravitational and Cosmological Theories 

In Subsec. 18.51 we showed that to construct a PT-symmetric model of quantum 
electrodynamics, one need only replace the electric charge e by ie and then replace 
the vector potential A^^ by an axial-vector potential. The result is a non-Hermitian but 
PT-symmetric Hamiltonian. An interesting classical aspect of this model is that the 
sign of the Coluomb force is reversed, so that like charges feel an attractive force and 
unlike charges feel a repulsive force. 

One can use the same idea to construct a PT-symmetric model of massless spin-2 
particles (gravitons). One simply replaces the gravitational coupling constant G by iG 
and then requires the two-component tensor field to behave like an axial tensor under 
parity reflection. The result is a non-Hermitian "PT-symmetric Hamiltonian, which at 
the classical level describes a repulsive gravitational force. It would be interesting to 
investigate the possible connections between such a model and the notion of dark energy 
and the recent observations that the expansion of the universe is accelerating. 

The connection discussed in Subsec. 12.61 between reflectionless potentials and VT 
symmetry may find application in quantum cosmology. There has been much attention 
given to anti-de Sitter cosmologies |133j and de Sitter cosmologies |134l 1135] . In the 
AdS description the universe propagates reflectionlessly in the presence of a wrong-sign 
potential {—x^, for example). In the dS case the usual Hermitian quantum mechanics 
must be abandoned and be replaced by a non-Hermitian one in which there are 'meta- 
observables'. The non-Hermitian inner product that is used in the dS case is based 
on the CVT theorem in the same way that the CVT inner product is used in PT- 
symmetric quantum theory [38]. The condition of reflectionless, which is equivalent 
to the requirement of PT symmetry, is what allows the wrong-sign potential to have 
a positive spectrum. Calculating the lowest energy level in this potential would be 
equivalent to determining the cosmological constant [136] . 
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8.8. Classical Field Theory 

The procedure for constructing PT-symmetric Hamiltonians is to begin with a 
Hamiltonian that is both Hermitian and VT symmetric, and then to introduce a 
parameter e that extends the Hamiltonian into the complex domain while maintaining 
its VT symmetry. This is the procedure that was used to construct the new kinds 
Hamiltonians in (|T2|) . We can follow the same procedure for classical nonlinear wave 
equations because many wave equations are VT symmetric. 

As an example, consider the Korteweg-de Vries (KdV) equation 



To demonstrate that this equation is VT symmetric, we define a classical parity 
reflection V to be the replacement x — > —a;, and since u = u{x,t) is a velocity, the 
sign of u also changes under V: u ^ —u. We define a classical time reversal T to be 
the replacement t — >• — t, and again, since m is a velocity, the sign of u also changes 
under T: u ^ —u. Following the quantum-mechanical formalism, we also require that 
i —y —i under time reversal. Note that the KdV equation is not symmetric under P or T 
separately, but that it is symmetric under combined VT reflection. The KdV equation is 
a special case of the Camassa-Holm equation |137j . which is also VT symmetric. Other 
nonlinear wave equations such as the generalized KdV equation Ut + u^u^ + u^xx = 0, 
the Sine-Gordon equation Uu — u^x + g sinu = 0, and the Boussinesq equation are VT 
symmetric as well. 

The observation that there are many nonlinear wave equations possessing VT 
symmetry suggests that one can generate rich and interesting families of new complex 
nonlinear PT-symmetric wave equations by following the same procedure that was used 
in quantum mechanics and one can try to discover which properties (conservation laws, 
solitons, integrability, stochastic behavior) of the original wave equations are preserved 
and which are lost. One possible procedure for generating new nonlinear equations from 
the KdV equation is to introduce the real parameter e as follows: 



Various members of this family of equations have been studied in Ref. [138] . Of course, 
there are other ways to extend the KdV equation into the complex domain while 
preserving VT symmetry; see, for example Ref. [139j . 

9. Final Remarks 

In this paper we have shown how to extend physical theories into the complex domain. 
The complex domain is huge compared with the real domain, and therefore there are 
many exciting new theories to explore. The obvious potential problem with extending 
a real theory into complex space is that one may lose some of the characteristics that 
a valid physical theory must possess. Thus, it is necessary that this complex extension 
be tightly constrained. We have shown that the essential physical axioms of a quantum 



Ut + UUx + M; 



0. 



(217) 
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theory are maintained if the complex extension is done in such a way as to preserve VT 
symmetry. 

The complex theories that we have constructed are often far more elaborate and 
diverse than theories that are restricted to the real domain. Upon entering the complex 
world we have found a gold mine of new physical theories that have strange and 
fascinating properties. We have just begun to study the vast new panorama that has 
opened up and we can hardly begin to guess what new kinds of phenomena have yet to 
be discovered. 

As an Ulam Scholar, CMB receives financial support from the Center for Nonlinear Studies at the Los 
Alamos National Laboratory and he is supported by a grant from the U.S. Department of Energy. 
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